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Abstract 

The effort to resolve data quality issues and improve on the initial data evaluation methodologies of the SeaWiFS 
Project was an extensive one. These evaluations have resulted, to date, in three major reprocessings of the entire 
data set where each reprocessing addressed the data quality issues that could be identified up to the time of 
the reprocessing. Three volumes of the SeaWiFS Postlaunch Technical Report Series (Volumes 9, 10, and 
11) are needed to document the improvements implemented since launch. Volume 10 continues the sequential 
presentation of postlaunch data analysis and algorithm descriptions begun in Volume 9. Chapter 1 of Volume 10 
describes an absorbing aerosol index, similar to that produced by the Total Ozone Mapping Spectrometer 
(TOMS) Project, which is used to flag pixels contaminated by absorbing aerosols, such as, dust and smoke. 
Chapter 2 discusses the algorithm being used to remove SeaWiFS out-of-band radiance from the water-leaving 
radiances Chapter 3 provides an itemization of all significant changes in the processing algorithms for each 
of the first three reprocessings. Chapter 4 shows the time series of global clear water and deep-water (depths 
greater than 1,000 m) bio-optical and atmospheric properties (normalized water-leaving radiances, chlorophyll, 
atmospheric optical depth, etc.) based on the eight-day composites as a check on the sensor calibration stability 
Chapter 5 examines the variation in the derived products with scan angle using high resolution data aroun 
Hawaii to test for residual scan modulation effects and atmospheric correction biases. Chapter 6 provides 
a methodology for evaluating the atmospheric correction algorithm and atmospheric derived products using 
ground-based observations. Similarly, Chapter 7 presents match-up comparisons of coincident satellite and in 
situ data to determine the accuracy of the water-leaving radiances, chlorophyll a, and (490) products. 


PROLOGUE 

The SeaWiFS Project Calibration and Validation Team 
(CVT) is responsible for the overall quality of the data 
products and for verifying the processing code. Volume 
38 of the SeaWiFS Technical Report Series (Prelaunch) 
outlined the prelaunch quality control strategy. Since Sea- 
WiFS began routine data processing in September 1997, 
the CVT has constantly worked to resolve data quality is- 
sues and improve on the initial data evaluation methodolo- 
gies. These evaluations resulted in three major reprocess- 
ings of the entire data set (February 1998, August 1998, 
and May 2000). Each reprocessing addressed the data 
quality issues that could be identified up to the time of 
each reprocessing. 

The number of chapters (21) needed to document this 
extensive work in the SeaWiFS Postlaunch Technical Re- 
port Series requires three volumes: Volumes 9, 10, and 
11. Volume 10 continues the sequential presentation of 
postlaunch data analysis and algorithm descriptions be- 
gun in Volume 9. The chapters describe the various data 
quality issues, analyses, and algorithm improvements that 
were developed through the third reprocessing. The data 
evaluations after the third reprocessing indicate that the 
data products, in most situations including Case-2 wa- 
ters, are within the prelaunch Case-1 water accuracy goals. 
Nonetheless, it is expected that other improvements and 
new geophysical data products will be developed in the 
future which will require additional reprocessings. The 
SeaWiFS Project Office will remain dedicated to provid- 
ing better products and to the documentation of future 
analysis and algorithm improvement studies. 


A short synopsis of each chapter in this volume is given 
below. 

1 . The Description of the SeaWiFS 

Absorbing Aerosol Index 

The method of detecting absorbing aerosols from Sea- 
WiFS is described in this chapter. This technique uses a 
term defined as the absorbing aerosol index which mea- 
sures the wavelength-dependent change in Rayleigh scat- 
tered reflectance due to aerosol absorption relative to a 
pure Rayleigh atmosphere. The sensitivities of the Sea- 
WiFS absorbing aerosol index to various types of absorb- 
ing and non-absorbing aerosols were shown using radiative 
transfer codes. Discussions of separating the signature of 
absorbing aerosols from that of phytoplankton on the spec- 
tral dependence of SeaWiFS reflectances are also included. 

2. Analyses of the SeaWiFS 
Spectral Band-Pass Effects 

This chapter describes an effort to study the effects of 
the SeaWiFS spectral band pass on the retrieved normal- 
ized water-leaving radiances and ocean near-surface chlo- 
rophyll concentrations. SeaWiFS routinely provides ocean 
near-surface optical and microphysical property data, how- 
ever, these retrievals include the effects of the relatively 
wide SeaWiFS spectral band passes. Because the SeaWiFS 
bio-optical algorithm was derived using in situ measure- 
ments, which were acquired at the SeaWiFS nominal band 
center wavelengths, it was necessary to assess the effects 
of the SeaWiFS spectral band pass on the retrieved chlo- 
rophyll concentrations. 
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3. Changes Made in the Operational 

SeaWiFS Processing 

This chapter discusses the major changes made in the 
SeaWiFS processing for the second reprocessing in August 
1998 and the third reprocessing in May 2000. Each major 
change in the processing programs and tables is presented 
with a description of the effects of the change. For the 
second reprocessing, a review of the major changes to the 
data is presented using an eight-day period of data as an 
example. The effect of the changes for the third reprocess- 
ing are also examined and compared to the results of the 
second reprocessing using two eight-day periods of data. 
The changes have increased the usefulness of the data in 
both oceanographic and atmospheric applications. 

4. SeaWiFS Global Clear-Water Analysis 

The SeaWiFS CVT made a comparison of global dear- 
water radiances retrieved by SeaWiFS with normalized 
water-leaving radiances measured by the Marine Optical 
Buoy (MOBY) as a check on the accuracy and stability of 
the vicarious calibration of SeaWiFS. The procedures and 
results of this comparison are described in this chapter. 

5. Along-Scan Effects in SeaWiFS Data 

The SeaWiFS CVT has looked for along-scan effects 
in SeaWiFS data by examining mean radiances of local 
area coverage (LAC) scenes over MOBY and by analyz- 
ing the ratio of SeaWiFS-to-MOBY (S:M) match-up data. 
Analyses of the SeaWiFS along-scan data and S:M ratios 
show decreases in the water- leaving radiances retrieved by 
SeaWiFS as the scan angles or optical paths of the ob- 
servations increase. These two analyses, which are inde- 
pendent of each other, both point to an overcorrection of 
the SeaWiFS data by the SeaWiFS atmospheric correction 
algorithm. 


6- SeaWiFS Aerosol Optical 

Thickness Match-up Analyses 

In this chapter, a match-up procedure is described 
which compares the retrieved SeaWiFS aerosol optical 
thicknesses with data from in situ measurements. The 
aerosol optical thickness at 865 nm is a by-product of the 
SeaWiFS atmospheric correction and is routinely retrieved 
from SeaWiFS measurements. This work is part of the 
SeaWiFS calibration and validation efforts in studying the 
aerosol optical properties over the ocean, thereby, validat- 
ing aerosol models used in the atmospheric correction of 
ocean color sensors. The aerosol model is an integral part 
of the SeaWiFS atmospheric correction. The SeaWiFS 
aerosol retrieval algorithm, the data acquisitions from both 
SeaWiFS and the in situ measurements, and the match-up 
procedure are described. Finally, some preliminary com- 
parison results are presented and discussed. 

7. Normalized Water- Leaving 

Radiance and Chlorophyll a 
Match-up Analyses 

Validation of SeaWiFS requires the use of in situ (field 
collected) data sets. The SeaWiFS and the Sensor Inter- 
comparison and Merger for Biological and Interdisciplinary 
Oceanic Studies (SIMBIOS) Projects have sponsored nu- 
merous Pis to collect in situ optics and chlorophyll data 
for the purpose of comparing values to those derived from 
the SeaWiFS instrument. The match-up design described 
here uses field data stored in SeaWiFS Bio-optical Archive 
and Storage System (SeaBASS), match-up analysis soft- 
ware, and a plotting and statistics package to validate the 
primary SeaWiFS derived ocean products, i.e., normalized 
water-leaving radiance, chlorophyll a concentration, and 
the diffuse attenuation coefficient. 
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Chapter 1 


The Description of the SeaWiFS 
Absorbing Aerosol Index 

N. Christina Hsu and Wayne D. Robinson 
SAIC General Sciences Corporation 
Beltsville , Maryland 

Sean W. Bailey 
Futuretech Corporation 
Greenbelt , Maryland 

P. Jeremy Werdell 

Science Systems and Applications , Incorporated 
Lanham, Maryland 

Abstract 

The method of detecting absorbing aerosols from SeaWiFS is described in this chapter. This technique uses 
a term defined as the absorbing aerosol index which measures the wavelength-dependent change in Rayleigh 
scattered reflectance due to aerosol absorption relative to a pure Rayleigh atmosphere. The sensitivities of e 
SeaWiFS absorbing aerosol index to various types of absorbing and non-absorbing aerosols were shown using 
radiative transfer codes. Discussions of separating the signature of absorbing aerosols from that of phytoplankton 
on the spectral dependence of SeaWiFS reflectances are also included. 


1.1 INTRODUCTION 

The presence of absorbing aerosols leads to errors in de- 
termining water-leaving radiances and, therefore, retriev- 
ing chlorophyll a concentrations. In this chapter, a method 
is described to distinguish absorbing aerosols from non- 
absorbing aerosols over the ocean, using an absorbing aero- 
sol index (AI) similar to the one developed for the Total 
Ozone Mapping Spectrometer (TOMS) sensor (Hsu et al. 
1996 and Herman et al. 1997). 


1.2 AEROSOL INDEX 

In the case of SeaWiFS, AI is determined from the 412 
and 490 nm remote sensing reflectances, R rs W, which have 
weak or negligible dependence on ozone and water vapor 
absorption. A compact notation for the reflectance ratio, 
, is used to represent R rs {K)/ Rrsi^j) i n the definition 
for'AI: 


I = -100 


l°gio (-^ 490 ) 


l°gl0 (^490 ) 


(1) 


atmosphere containing only large nonabsorbing particles 
and Rayleigh scattering molecules. Essentially, I is a mea- 
sure of the wavelength-dependent change in Rayleigh scat- 
tered reflectance due to aerosol absorption relative to a 
pure Rayleigh atmosphere. 

AI is defined so that positive values generally corre- 
spond to absorbing aerosols and negative values to non- 
absorbing aerosols (Hsu et al. 1996 and Herman et al. 
1997). The AI measurements are linearly proportional 
to the aerosol optical thickness (AOT) within the range 
of AOT observed from ground-based sun photometer in- 
struments in regions covered by burning biomass smoke or 
wind-blown African dust (Hsu et al. 1999). 

In the following sections, the sensitivity of the AI to 
the detection of various types of absorbing aerosols over 
blue water is discussed. The effects of phytoplankton on 
absorbing aerosol detection and how to account for them 
is also described. Finally, the potential applications of the 
SeaWiFS AI to general algorithmic issues is outlined. 

1.3 SENSITIVITY STUDIES 


where / is the aerosol index, R 490 * s ^he top of the at- 
mosphere (TOA) reflectance ratio measured by SeaWiFS, 
and R 490 is the reflectance ratio calculated assuming an 


Sensitivity studies were performed for various types of 
absorbing and nonabsorbing aerosols using radiative trans- 
fer codes (Dave 1972). To use aerosol properties represen- 
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Table 1 . The characteristics of the aerosol properties used in the model simulations for various types of 
aerosols. The real and imaginary parts of the refractive index are shown in complex notation with i denoting 
the square root of -1 (v-1). 


Aerosol 

Model 

r oM 

a' 

Refractive Index 
(412 nm) 

Refractive Index 
(490 nm) 


Strongly Absorbing Dust 

0.25 

2.20 

1.57-0.015i 

1.57 — O.OIOz 


Less Absorbing Dust 

0.12 

2.20 

1.57-0.015* 

1.57-0.010i 


Sea Salt 

2.00 

2.15 

1.38-0.000* 

1.38-0.000* 


Darker Smoke 

0.14 

1.45 

1.55-0.040i 

1.55-0.035* 


Whiter Smoke 

0.14 

1.45 

1.55—0.023* 

1.55-0.020* 


Sulfate 

0.07 

2.03 

1.43—0.000* 

1.43-0.000* 




^*8* The relationship between the SeaWiFS AI and AOT at 490 nm for various types of tropospheric 
aerosols over blue waters. The corresponding single scattering albedo (u/ a ) is shown next to the aerosol 
model. The height of the aerosol layer is assumed to be at 3 km. The viewing geometry is for the condition 
of the solar zenith angle of 40°, satellite zenith angle of 52°, and relative azimuth angle of 120°. 
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tative of smoke from biomass burning and dust particles 


from sandstorms, values based on the results of in situ 
aerosol measurements (Tegen and Lacis 1996, and Remer 
et al. 1996) were employed in the model simulations. A 
log-normal particle size distribution was assumed in this 
study. The mode radius, r 0 (/x)> and width, <r', of the par- 
ticle size distribution and the real and imaginary part of 
the refractive index at 412 and 490 nm for the aerosol mod- 
els used in these calculations are listed in Table 1. 

Figure 1 illustrates the relationship between the Sea- 
WiFS AI and AOT over blue waters for the aerosol mod- 
els listed in Table 1. The AOT corresponds to values 
at 490 nm. The mean aerosol altitude is assumed to be 
3.0 km. The results in this figure are calculated for a solar 
zenith angle of 40°, a satellite zenith angle of 52°, and a 
relative azimuth angle of 120°. It is apparent that the slope 
of the AI versus optical depth changes with aerosol type 
(i.e., single scattering albedo, u; a ). Whiter smoke (e.g., 
a/ 0 = 0.91) can only be detected at larger aerosol optical 
thicknesses. On the other hand, more strongly absorbing 
dust aerosols (e.g., t o a = 0.80) can be detected at smaller 
optical thicknesses using this absorbing aerosol index tech- 
nique. 

1.4 CASE-2 WATER CONSIDERATIONS 

The signals over an ocean surface due to highly absorb- 
ing phytoplankton can mimic the spectral characteristics 
of absorbing aerosols in the remotely sensed reflectance. 
To help distinguish signals due to absorbing aerosols from 
those due to ocean surfaces with high chlorophyll content, 
the SeaWiFS 765 nm channel is used. The reflectance mea- 
sured at this channel is greatly influenced by gas absorp- 
tion associated with the oxygen A-band. 1 he magnitude of 
such oxygen absorption depends on the effective altitude of 
the reflecting surface (e.g., ocean surface, cloud, or aerosol 
layer surface). It is assumed that the general sources of ab- 
sorbing aerosols are from land, and the height of a result- 
ing aerosol layer is not close to ocean surfaces. By selecting 


Schmidt, Wang, Bailey, and Werdell 

a threshold on the magnitude of the oxygen A-band sig- 
nal, low-altitude aerosols (i.e., sea-salt) can be screened 
out, and conditions where a sufficient amount of absorb- 
ing aerosols are present to retrieve their properties can be 
determined. 

When there is a intermediate or high amount of aerosol 
loading, the detection of absorbing aerosols is less sensitive 
to surface properties. For a thin aerosol layer over Case-2 
water, however, it is very difficult to differentiate between 
absorbing and nonabsorbing aerosols. There is also a sig- 
nificant uncertainty in modeling the reflectance of various 
types of constituents over the ocean at 412 and 490 nm. 
Instead of trying, therefore, to model the effects, a chlo- 
rophyll content-based threshold is used as a discriminator. 
A tighter constraint on this threshold was employed when 
determining absorbing aerosol contamination for Case-2 
water than for the Case-1 water. To determine the thresh- 
old level, a simple chlorophyll index was calculated using 
the SeaWiFS blue and green channels to provide a rough 
estimate on the degree of deviation of target water from 
blue water. When a higher level of chlorophyll content is 
indicated, a larger value of AI is required to pass the con- 
fidence level for detecting absorbing aerosols. As a result, 
the detection limit for absorbing aerosols is higher over 
Case-2 water than that over Case-1 water. 

The combination of procedures mentioned above has 
provided reasonably good results to distinguish absorbing 
aerosols from the phytoplankton signature using the re- 
motely sensed reflectance. 

1.5 POTENTIAL APPLICATIONS 

The SeaWiFS AI provides a fast way to identify the ab- 
sorbing aerosols on a pixel- by-pixel basis. It can be used 
in the current operational algorithm as a flag to filter out 
chlorophyll and water-leaving radiance retrievals contam- 
inated by absorbing aerosols. It will also aid the aerosol 
model selection process once the absorbing aerosol models 
are included in the retrieval algorithm. 
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Chapter 2 


Analyses of the SeaWiFS 
Spectral Band-Pass Effects 

Menghua Wang 

University of Maryland, Baltimore County 
Baltimore, Maryland 

Bryan A. Franz and Robert A. Barnes 
SAIC General Sciences Corporation 
Beltsville , Maryland 

Abstract 

This chapter describes an effort to study the effects of the SeaWiFS spectral band pass on the retrieved nor- 
malized water-leaving radiances and ocean near-surface chlorophyll concentrations. SeaWiFS routinely provides 
ocean near-surface optical and microphysical property data, however, these retrievals include the effects of the 
relatively wide SeaWiFS spectral band passes. Because the SeaWiFS bio-optical algorithm was derived using 
in situ measurements, which were acquired at the SeaWiFS nominal band center wavelengths, it was necessary 
to assess the effects of the SeaWiFS spectral band pass on the retrieved chlorophyll concentrations. 


2.1 INTRODUCTION 


spectral response functions. By defining 


SeaWiFS, which was successfully launched on 1 Au- 
gust 1997, is one of only a few satellite instruments having 
complete prelaunch band spectral response measurements 
covering wavelengths from 380-1,150 nm for the eight Sea- 
WiFS bands. One of the stated goals of the mission is to 
estimate water-leaving radiances to an accuracy of ±5% 
(Hooker et al. 1992). Figures 2a-d provide the spectral 
response function (SRF) for SeaWiFS bands 2, 3, 5, and 
6. SeaWiFS not only has in-band response structures, but 
it also has significant sensor out-of-band contributions. 

In ocean color remote sensing, the sensor- measured ra- 
diance at the top of the ocean-atmosphere system, mea- 
sured at a wavelength A, can be written as, 

L t ( A) = L r ( A) +L a (A) + L ra (A) + t(X)L w (X), (2) 

where L r (A), L a (A), and L ra (A) are the radiance con- 
tributions from the multiple scattering of air molecules, 
aerosols, and Rayleigh- aerosol interactions, respectively; 
L\v{ A) is the water- leaving radiance, which is the desired 
quantity for relating the oceanic near-surface physical and 
bio-optical properties; and t( A) is the atmospheric diffuse 
transmittance which accounts for the effects of propagat- 
ing Lw{ A) from the sea surface to the TOA (Wang 1999a, 
and Yang and Gordon 1997). 

Because of sensor spectral band responses, the average 
radiance measured by SeaWiFS is weighted by the sensor 


/ m , V = lL(X)S i (\)d\ 
(L(A,) > JS t (X)dX ’ 


(3) 


where S*(A) is the SeaWiFS spectral response function for 
band i at a nominal center wavelength A l} (3) can then be 
rewritten as 


(Lt(X t )) — (L r (A,)) -f- (Z/ a (A t ) + L ra (A t )} 

+ (t(Xi) Lw(Ai)). (4) 

In (4), i = 1-8 selects the SeaWiFS eight spectral bands 
and (L t ( A,)} is the radiance at the TOA measured by the 
SeaWiFS instrument. The Rayleigh contribution (L r (A t )), 
as well as the aerosol and Rayleigh-aerosol interaction term 
{L a ( A t ) -F L Ta {X % )), can be estimated using methods out- 
lined in Gordon (1995). The Gordon and Wang atmo- 
spheric correction algorithm can then be applied (Gordon 
and Wang 1994a). 

Gordon (1995) showed the averaged Rayleigh and ozone 
optical thicknesses, as weighted by the solar irradiance and 
the sensor band SRF in computing the atmospheric diffuse 
transmittance, can be used to rewrite the last term in (4): 

(t(X t ) L W (X,)) X t(Xi)t 0 (X % ) cos 6 0 {Lwn(X 1 )), (5) 

where Lwn{ At) is the normalized water-leaving radiance 
for a given band, and fo(A,) and £(A*) are the atmospheric 
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Fig. 2. (cont.) The SeaWiFS spectral response functions for SeaWiFS: d) band 6 (670 nm). 



Fig. 3. The spectral distribution of the normalized water-leaving radiance Lwn( A) for a typical Case-1 water 
for various pigment concentrations, derived from Gordon et al. (1988). 


diffuse transmittance at the solar and instrument view- 
ing direction, respectively. SeaWiFS, therefore, routinely 
reports the weighted normalized water-leaving radiances, 
and it is the ratios of these weighted values that are used 
to compute the oceaninc near-surface chlorophyll concen- 
trations (O’Reilly et al. 1998). 

2.2 SPECTRAL BAND-PASS EFFECTS 

To assess the effects of SeaWiFS spectral band respon- 
ses on the derived normalized water-leaving radiance, it 
is necessary to derive a consistent set of results for both 
the nominal center wavelengths, Lwn(K), and the SRF 
weighted, (Xvw(A,)), cases. The spectral distribution of 


LwnW was, therefore required at various chlorophyll con- 
centrations. Figure 3 provides typical L W nW (Case-1 wa- 
ter) spectral distributions obtained from bio-optical model 
calculations (Gordon et al. 1988) for four pigment concen- 
tration values, C a +Cp = 0.01, 0.1, 0.5, and l.Omgm' 3 , 
where C a +C P is the sum of chlorophyll a and phaeophytin a 
concentrations, respectively. 

The Fig. 3 values were used to compute (L W n( A*)) 
and to estimate the differences between (Tw 7 v(Ai)) and 
Lwn( Ai). It was found that the values in SeaWiFS bands 
4 and 5 (green bands) for clear ocean water (low pigment 
concentrations) overestimated the Lwn{K) values by ap- 
proximately 5% (Fig. 4). Because the SeaWiFS band 5 
(555 nm) value was used in deriving the pigment concen- 
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450 500 550 600 650 

Wavelength [nm] 


Fie. 4. The difference in percent between values of (Lwjv(A»)) and L WN (\ t ) at the SeaWiFS spectra an s 
for pigment concentration in milligrams per cubic meter of a) 0.01, b) 0.1, and c) 0.5. Each figure also gives 
the values of (lm(A|)} (scale noted in the right side of figure). 
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Wavelength [nm] 


Fig. 4. (cont.) The difference in percent between values of (Lvw(Ai)} and Lwn{\) at the SeaWiFS spectral 
bands for chlorophyll concentration in milligrams per cubic meter of d) 1.0. The figure also gives the values 
of (Lwn{ At)) (scale noted in the right side of figure). 


trations, it would be useful if the SeaWiFS SRF effects 
could be removed and converted to the L W n{ K) value. 

2.3 CORRECTIONS 

To convert (L W n{K)) to L WN {\ % ) at the SeaWiFS 
spectral bands, the spectral distribution of Lwn{K) is 
needed; this depends on the ocean near-surface optical 
and microphysical properties. The spectral distribution of 
Lwn(K)> however, can not be predicted a priori . Because 
SeaWiFS retrieves (L^ivfA,)) at the six spectral bands 
(412, 443, 490, 510, 555, and 670 nm), a first approxi- 
mation for Lwat(A) can be derived by interpolating the 
(Lwn(K)) values, and then using this measured spectral 
distribution to derive correction coefficients. 

Using this approach, an iterative procedure can be de- 
fined to convert the SeaWiFS (L W n{K)) values to band- 
centered wavelength values, L' WN ( X t ): 

1. SeaWiFS retrieved (Lwn{ A,)) can be approximated 
as Lwn{ Ai), i.e., Lwn{ A*) « (Lwn{ A t )}. 

2. Using Lwn{ A,) values, the average of the parameter 
(Z/( A,)) as weighted by the sensor spectral response 
function, |Lw/v(Ai)), can be derived using (3). The 
ratio values, r(A t ) (i.e., the ratio value of Lwn at 
the SeaWiFS band-centered wavelength versus the 
average values weighted by the SeaWiFS spectral 
response function), can be estimated with 


r <\ \ = 

1 %) (Lwn(K))' 


(6) 


3. Finally, the normalized water-leaving radiance Sea- 
WiFS band-centered wavelength, U WN (Xi)^ can be 
obtained with 

U WN (\ X ) = r(\i) (Lwn(K))- (7) 

Going from step 2 to step 3 is an iterative process. In 
testing this approach, it was found that one iteration was 
usually accurate. It should be noted that to compute the 
required integrals, a log-linear interpolation was used to 
convert the six SeaWiFS {L W n{ A*)) values to a relatively 
continuous wavelength distribution covering the SeaWiFS 
spectral range from 380-1,150 nm. 

Figures 4a-d provide results of the correction scheme 
for pigment concentrations of 0.01, 0.1, 0.5, and 1.0 mg 
m~ 3 . Although the effects are negligible in the SeaWiFS 
blue bands, the corrections were performed for all six visi- 
ble bands to test the efficacy of the approach. In generating 
Figs. 4a-d, the spectral distribution of Lwn{ A), as shown 
in Fig. 3, was used. Note that in Figs. 4a-d, there are 
three curves in each figure: 

a) The SeaWiFS (Lwjv(A,)) values (the scale is noted 
on the right side of the figure); 

b) The difference between and Lwn{K) 

without corrections; and 

c) The same as in item b, but with corrections. 

The correction scheme works quite well. Differences be- 
tween true and retrieved values of Z/wvv(Ai) are reduced, 
in particular, in the SeaWiFS green and red bands. The 
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correction scheme was also applied to in situ data sets ac- 
quired from MOBY in Hawaii, and similar results were 
obtained. 

This correction scheme was implemented into the chlo- 
rophyll match-up analysis, where the SeaWiFS derived 
chlorophyll values were compared with co-located in situ 
measurements. The results were plotted as SeaWiFS re- 
trieved values versus the in situ measurements, and a lin- 
ear least-squares fit was then obtained. Table 2 gives the 
values of intercept (fr), slope (m), and the correlation coef- 
ficient ( R ) before and after the SeaWiFS band-pass effects 


corrections. Note that the perfect match corresponds to 
values of b = 0, m = 1, and R = 1. 

Table 2. Values of coefficients in the linear least- 
squares fit for the chlorophyll match-up analyses 
(SeaWiFS versus in situ measurements) before and 
after the SeaWiFS band-pass effects corrections. 


Case 

b 

m 

R 

Before 

After 

0.1093 

0.0824 

1.1272 

1.0266 

0.8665 

0.8504 
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Abstract 

This chapter discusses the major changes made in the SeaWiFS processing for the second reprocessing in 
August 1998 and the third reprocessing in May 2000. Each major change in the processing programs and tables 
is presented with a description of the effects of the change. For the second reprocessing, a review of the major 
changes to the data are presented using an eight-day period of data as an example. The effect of the changes 
for the third reprocessing are also examined and compared to the results of the second reprocessing using two 
eight-day periods of data. The changes have increased the usefulness of the data in both oceanographic and 
atmospheric applications. 


3.1 INTRODUCTION 

This chapter discusses the major changes made in the 
SeaWiFS level-1, -2, and -3 processing which led to the 
third reprocessing. Ever since the first SeaWiFS reprocess- 
ing of 2 January 1998, SeaWiFS Project personnel have 
been evaluating the results of the processing and have dis- 
covered many problems and improvements. Some prob- 
lems were fixed in the operational software, but many im- 
provements which required recalibration could not be ap- 
plied without initiating a complete reprocessing. These 
changes were accumulated for a period of approximately 
six months before applying them to all the data in the 
second reprocessing of the data set. During this period, 
significant changes in the calibration of SeaWiFS bands 
7 and 8 (Eplee and McClain 2000a) made it imperative 
to include the calibration and other changes into the op- 
erational processing, so the quality of the data could be 
maintained. 

The second reprocessing started on 14 August 1998. 
At that time, the new software was used to reprocess all 
of the data collected from the start of the mission and 
to process the real-time data as it was acquired. Since the 
second reprocessing, a number of minor enhancements and 


repairs were made to the operational processing software. 

Several remaining problems, which were recognized af- 
ter the second reprocessing, were extensively investigated; 
many improvements were developed to make the resulting 
data more useful and to broaden the product suite. All 
these changes were incorporated in the third reprocessing, 
which occurred in May 2000. 

The next section (Sect. 3.2) discusses the major changes 
made in the software and tables for the second reprocess- 
ing, and the impact of those changes on the products. A 
brief summary of the combined effect of those changes is 
discussed in Sect. 3.2.4 and the major changes included 
in the third reprocessing are presented in Sect. 3.3. Sec- 
tion 3.3.6 presents an analysis of the changes in the derived 
products from the second to the third reprocessings, which 
is followed by a summary of the changes and a discussion 
of future improvements (Sect. 3.4). 

3.2 THE SECOND REPROCESSING 

In the approximately six months between the first and 
second SeaWiFS reprocessings, many changes were made 
in the level -2 and -3 processing software. The most signif- 
icant changes in that period are shown in Table 3 and are 


12 


McClain, Barnes, Eplee, Jr., Franz, Hsu, Patt, Pietras, Robinson, Schieber, Schmidt, Wang, Bailey, and Werdell 


T^ble 3. Major changes and their observable effects in the second reprocessing. 

Change 

Observable Effect 

Calibration: 

Prelaunch Calibration Change 
Time dependence of bands 7 and 8 
Revised vicarious calibration 

None 

Better temporal stability 
None 

Atmospheric Correction: 

New transmittance tables 
Wider range of aerosol models 

Removed Lwn as / (scan angle) 

Reduced number of points outside model range 

Other: 

Updated coefficients for the chlorophyll a 
algorithm 

Added test for high reflectance ratios 
Tests to decrease chlorophyll speckles 
Improved Lwn computation 
Noise resistant dark count 
Improved solar radiance computation 
More accurate anchor point grid 
Reduced whitecap radiance 
Binned higher chlorophyll 
Wrote out actual t values 

Increased low chlorophyll a concentrations and 
decreased high chlorophyll a concentrations 
Removed spurious chlorophyll 
Some reduction on speckles 
Less than 1% change in Lwn 
No striping in the data 
None 
None 

Reduced data loss 

Retained larger chlorophyll range 

Better diagnostics on e values 


discussed briefly here. The details of changes related to 
the calibration, atmospheric correction, and the new chlo- 
rophyll a algorithm are discussed elsewhere in this three- 
volume set, Volumes 9-11, of the SeaWiFS Postlaunch 
Technical Report Series . 

3.2*1 Calibration Related Changes 

In the period leading up to the second reprocessing, the 
calibration of the SeaWiFS instrument became more rig- 
orous, with a more detailed band 7 calibration and better 
detector characterization. Three changes associated with 
the calibration are presented here. 

3. 2. 1.1 Prelaunch Calibration Changes 

The first reprocessing was conducted using a prelimi- 
nary set of SeaWiFS calibration values, which were deter- 
mined in January and April 1997. When an analysis of 
the calibration data was completed (Johnson et al. 1999), 
the calibration was modified to include these updates, and 
was used for the second reprocessing. The changes in the 
total radiances were less than 0.5%. 

3.2. 1.2 Time Dependent Calibration 

In 1998, the time series of solar and lunar calibration 
measurements, which began in September 1997, showed 
that the radiometric responses of bands 7 and 8 were chang- 
ing as a function of time. The greatest change occurred in 
band 8, amounting to a 4-5% decrease in the total radi- 
ances. The calibration table used during the second re- 
processing was adjusted to compensate for these changes 


and was updated as needed thereafter in the operational 
processing. A detailed discussion of the calibration is pre- 
sented in Eplee and Barnes (2000) and Eplee et al. (2000). 

3.2. 1.3 Revised Vicarious Calibration 

Improvements were made in the vicarious calibration 
of band 7 relative to band 8, and for bands 1-6. Details 
of the vicarious calibrations are discussed in Robinson and 
Wang (2000) and Eplee and McClain (2000a), but a short 
summary is presented here. 

The band 7 vicarious calibration was accomplished by 
adjusting the band 7 gain factor until the ratio of the 
single-scattering aerosol radiances (the e value) in bands 
7 and 8, matched the expected ratio over the open ocean. 
This process was applied to a limited number of sites in the 
first reprocessing. The current band 7 calibration extends 
and improves this procedure by selecting open-ocean sites 
near Hawaii, evaluating the gain using the time series of 
SeaWiFS LAC data supporting the MOBY site, and using 
strict screening procedures for data quality. 

The calibration for bands 1-6 was derived by adjusting 
the gain factors for each band until the water-leaving radi- 
ances in those bands matched those measured by MOBY. 
In the second reprocessing, the calibration was improved 
with more MOBY observations. 

3.2.2 Atmospheric Correction Changes 

Wang (2000) presents a more detailed description of the 
changes made in the atmospheric correction. Two major 
changes are briefly described here. 
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3.2.2. 1 New Transmittance Tables 

The transmittance tables were corrected to include the 
Fresnel reflectivity of the air-water interface. The intro- 
duction of this correction removed a small scan-angle de- 
pendence in Lwn ■ The changes reduced Lwn in bands 
1-5 by 2% at the center of the scan. At the edge of the 
SeaWiFS global area coverage (GAC) scan, the Lwn was 
reduced by about 5%, and at the edge of a LAC or high res- 
olution picture transmission (HRPT) scan, the reduction 
was about 10%. 

3. 2. 2. 2 Wider Range of Aerosol Models 

In the process of performing the atmospheric correc- 
tion, the level -2 processing determines the atmospheric 
model most closely matching the t value, or ratio of aerosol 
radiances in bands 7 and 8. For about 40% of the points 
processed in the first reprocessing, the e value was below 
the lowest e value of the available models. In these cases, 
the atmospheric correction was forced to use the model 
with the lowest e value to derive the aerosol radiances for 
the other bands. Although this model was the best avail- 
able, it was not as correct as a model that would fit the 
observed e value. 

In the second reprocessing, the model suite was ex- 
panded to include oceanic aerosol models which contain 
lower e values. The new atmospheric model suite reduced 
the occurrence of points below the model range from 40% 
to about 20%. 

3.2.3 Other Changes 

Other changes made to the SeaWiFS processing affect 
both the level -2 and -3 algorithms. The changes correct 
minor problems, improve algorithms, or reduce noise in the 
products. 

3.2.3. 1 Improved Chlorophyll a Algorithm 

The relationship between Lwn and chlorophyll a was 
updated with additional in situ observations (Maritorena 
and O’Reilly 2000), which resulted in a new set of coeffi- 
cients for the equation relating the ratio of Lwn to chloro- 
phyll in the 490 and 555 nm bands. The new relationship 
increased chlorophyll values that are less than 0.03 mg m~ 3 
and decreased chlorophyll values that are greater than 1.0 
mg m“ 3 . 

It was found that under conditions of very large re- 
flectance ratios, (values greater than 10), the chlorophyll 
algorithm could return reasonable chlorophyll values. The 
chlorophyll failure flag was set in these instances, so these 
values would not be binned. A similar treatment was ap- 
plied in the third reprocessing. 

3.2. 3. 2 Decreased Chlorophyll a “Speckling” 

The SeaWiFS level -2 and -3 binned data contained 
isolated pixels with relatively high chlorophyll a values in 


fields with otherwise low chlorophyll values. Frequently, 
these so-called speckles were near cloud edges not masked 
as “cloud” pixels. Methods for removing these points (such 
as, decreasing the band 8 albedo threshold for cloud mask- 
ing, or expanding the stray light distance), caused a large 
number of good data points to be masked (and, thus, re- 
moved) in the process of removing a small number of bad 
points. 

An alternative method was used that masked many 
occurrences of speckles without masking (and removing) 
good data. A mask was applied to the pixels in which 
any of the bands 1-8 had zero or negative values after the 
radiances were corrected to remove the effects of ozone ab- 
sorption, whitecap radiance, and Rayleigh radiance. The 
masked data points were assigned the atmospheric correc- 
tion failure flag (which is usually a mask). The new test re- 
duced the occurrences of speckles, but did not completely 
eliminate them. This test also revealed cases where the 
whitecap correction to the radiances was too large (Sect. 
3.2.3. 7). 

3. 2. 3. 3 Improved Lyyjy Computation 

The Lw normalization step of the level -2 processing 
was improved to include a more realistic atmospheric at- 
tenuation than the attenuation used in the first reprocess- 
ing. In the first reprocessing, the normalization only used 
the Rayleigh scattering and ozone absorption to correct for 
the attenuation in the path from the sun to the surface. 
For the second reprocessing, the effect of aerosol scatter- 
ing was included, which slightly reduced the Lwn values 
in bands 1-5 by less than 1%. 

3. 2. 3. 4 Noise Resistant Dark Count Calculation 

Striping was noticed in the radiance data, the cause 
for which was traced to the dark count calculation used in 
the calibration of the data. The dark count is the digital 
count value measured by the detectors when looking at a 
dark surface and, thus, represents a zero radiance value in 
the calibration equation. The dark count is collected for 
every scan line of data, has a nominal value of 20 counts, 
and remains virtually constant throughout an entire data 
pass. In some of the data, the dark count had a tendency 
to change by one digital count from one line to the next, 
which is commonly referred to as jitter. 

The jitter in the dark counts caused the line-to-line 
striping in the radiances. To remove the jitter, the dark 
counts were averaged to estimate a dark count offset for 
the entire pass. With the onset of data encryption, a num- 
ber of dark count values that were very different with re- 
spect to the standard values from undecrypted lines, were 
included in the algorithm that computes the average dark 
count value. The resulting dark count average occasionally 
produced badly calibrated radiances and caused the entire 
scene to have radically different normalized water-leaving 
radiances than expected. 
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The dark count calculation was improved using the fol- 
lowing two steps: 

1) Remove from consideration any dark count values 
that are less than 5 counts or greater than 35 counts, 
and 

2) Take the median value as the dark count for the 
pass. 

The improvements were especially helpful for HRPT data 
processing, which frequently has trouble with noise at the 
start and end of a pass. On occasion, it also helped improve 
the quality of GAC and LAC data. 

3. 2 .3. 5 Improved Solar Radiance Computation 

A 0.1% error in the value of the correction in the solar 
irradiance as a function of the time of year, was discovered 
in the level -2 processing. The correction was improved to 
reduce this error to 0.001%. A 0.1% error in the solar ir- 
radiance roughly translates to approximately a 1.0% error 
in the Lwn values. 

3. 2. 3. 6 More Accurate Anchor Point Grid 

To speed up the calculation of processing parameters, 
such as, the ancillary data values and Rayleigh radiances, 
the level -2 processing program used a set of anchor points 
at which these quantities were computed exactly. The val- 
ues of the ancillary data at other points were computed 
by interpolating the anchor point data to that location. 
For the first reprocessing, the anchor point separation was 
eight pixels for LAC and HRPT data, and two pixels for 
GAC data. 

In some cases, the difference between the interpolated 
Rayleigh values with respect to the exact values was sig- 
nificant, e.g., as great as 0.015 mW cm -2 /im" 1 sr 1 in the 
Rayleigh radiance at 443 nm, which roughly translates to a 
1% uncertainty in the Lwn- O n average, this uncertainty 
was less than 0.002 mW cm' 2 /im" 1 sr" 1 . This uncertainty 
was reduced to less than 0.004 mW cm 2 /im 1 sr 1 in the 
second reprocessing by reducing the anchor point separa- 
tion by a factor of 2 (one pixel in the GAC data and four 
pixels in the LAC and HRPT data) to improve the accu- 
racy of the Rayleigh and other ancillary data used. (Note 
that in the third reprocessing, all calculations are done at 
every point). 

3.2.3. 7 Reduced Whitecap Radiance Effect 

Several instances occurred where a region of ocean was 
cloud-free, but it was masked out as an atmospheric cor- 
rection failure, because the ozone, whitecap, and Rayleigh- 
corrected radiances in some bands were negative (a mask 
condition intended to reduce speckles in the chlorophyll 
concentration). It was discovered that in very clear re- 
gions with high wind speeds (around 15ms" 1 ), the white- 
cap correction was so large in bands 7 and 8, the remain- 
ing radiance was smaller than the Rayleigh radiance which 
caused the region to be masked. 


A 75% reduction in the whitecap correction was im- 
plemented to solve this problem. In discussions with H. 
Gordon, this reduction factor seemed reasonable. In ad- 
dition, studies of the whitecap contributions (Gordon and 
Wang 1994b) showed a large amount of variability in the 
whitecap radiances under high wind conditions. The main 
effect of this change in the whitecap correction radiances 
was to allow previously overcorrected and masked areas to 
be processed. 

3. 2.3. 8 Binning of Higher Chlorophyll Values 

It was noticed that in the process of binning the level -2 
data, a large amount of high chlorophyll data was being 
excluded. The primary reason for this was that many pix- 
els had a valid chlorophyll value (i.e., the Lwn in bands 
3 and 5 were positive), but the Lwn values in bands 1 
or 2 were negative. In the first reprocessing, the presence 
of negative Lwn values in a measurement was one of the 
exclusion conditions for the binning. 

To bin more of these excluded chlorophyll values, a 
number of changes were made in the level -2 and -3 pro- 
grams. First, the level -2 program was modified to output 
the actual negative water-leaving radiances that were de- 
rived instead of a —1.0 value. This allows an investigator 
to get a better idea of how large the negative Lwn values 
actually are. The binning routine was modified to bin pix- 
els that have negative Lwn values in them, but to bin them 
as zeros. In addition, the chlorophyll algorithm failure flag 
(CHLOftl) was added as a binning flag, so points which had 
a chlorophyll algorithm failure (band 3 or 5 where Lwn 
less than zero or where chlorophyll values are greater than 
64) would not get binned. This set of changes restored 
much of the previously excluded high chlorophyll data. 

3. 2. 3. 9 Actual Epsilon Values 

In the first reprocessing, the aerosol correction algo- 
rithm in the level -2 processing code reported the e value 
for the nearest aerosol model in cases where the derived e 
value was actually outside the model range. This caused 
some uncertainty in the actual value of e, and it also tainted 
some calculations that used the average value of e. 

The aerosol correction algorithm was modified to report 
and store the actual computed e value. This e value was 
also used in the aerosol calculations. 

3.2.4 Second Reprocessing Analysis 

In many cases, changes in Lwn values brought about 
by level -2 processing algorithm changes were compensated 
for when the vicarious calibration was applied. The only 
way to understand the effects of all the changes was to 
consider them after the calibration was implemented. The 
cumulative effect of all the above changes is discussed in 
this section. The effect of all the changes to the processing 
is studied using a sample of typical SeaWiFS data. 
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Chlorophyll 

Fig. 5. The SeaWiFS level-3 global chlorophyll a frequency distribution for an eight-day period (14-21 March 
1998) using the first and second reprocessing algorithms (dotted and solid curves, respectively). 


3.2. ^.1 Chlorophylls Changes 

Figure 5 shows the global frequency distribution of chlo- 
rophyll a resulting from the first and second reprocessings 
for an eight-day period. The second reprocessing chlo- 
rophyll values are higher in open ocean regions, such as, 
the South Pacific and Atlantic. Large regions that previ- 
ously had chlorophyll concentrations at the lowest value of 
0.01 mg m -3 with the first reprocessing, had their values 
increase to around 0.03 mgm -3 with the second reprocess- 
ing. This was mostly the result of the new chlorophyll 
algorithm, although some contribution resulted from the 
new transmittance tables and calibration. 

The increased chlorophyll values observed in the second 
reprocessing agreed better with observations of the lowest 
chlorophyll concentrations, i.e. , values of 0.01 mgm -3 are 
rarely observed. The distribution of chlorophyll was raised 
significantly in the second reprocessing. The peak value of 
chlorophyll increased from 0.02-0.05 mgm -3 and far fewer 
points had the lowest chlorophyll value of 0.01 mgm -3 . A 
number of observations still appeared below the lower limit 
of 0.01 mgm -3 in the second reprocessing. 

Another aspect of the second reprocessing was that the 
total number of filled bins in the time binned product in- 
creased by 5%. Many of the new points were in high chlo- 
rophyll areas — a direct result of the new binning strategy 
which retained a greater number of high chlorophyll val- 
ues (Sect. 3. 2. 3. 8). Many of the new points were in coastal 
regions and in the Baltic Sea. 

The second reprocessing significantly reduced the mag- 
nitude of the highest chlorophyll values (another direct 


effect of the new chlorophyll a algorithm). Chlorophyll 
concentrations near many coastal regimes, for example, 
showed decreases up to 30%. The second reprocessing 
changes also resulted in an increase of the Lwn values 
in bands 3 or 5, so that chlorophyll values could be deter- 
mined in more coastal areas. 

Both the increase in very low chlorophyll a values and 
the decrease in high chlorophyll a values were improve- 
ments that occurred in the new processing. 

3. 2.4* 2 Water- Leaving Radiance Changes 

An analysis of the global distribution of Lwn showed 
a general increase in the Lwn values with the second re- 
processing. For instance, the mean Lwn at 555 nm was 
increased by 0.06 from 0.26-0.32. The new mode of 0.27 
(versus 0.21 in the first reprocessing) had better agree- 
ment with the MOBY-derived clear-water Lwn value of 
0.254 (Eplee and McClain 2000b). 

In the blue bands, considerably less change was ob- 
served in the Lwn distribution. The mean Lwn at 412 nm 
increased only slightly from 1.53-1.56. In some isolated re- 
gions along the coast, Lwn actually decreased somewhat. 
The decreases may have been due to the inclusion of zero 
values of Lwn in the binning algorithm (Sect. 3.2. 3. 8). It 
was hoped that the second reprocessing changes would re- 
duce the number of negative Lwn values in bands 1 and 2. 
Although Lwn increased overall, there were still problems. 
The next section looks at the distribution of negative Lwn 
values in more detail. 
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Fig. 6. Percent occurrence of negative Lwn values at 412 nm for the eight-day period from 12-19 July 
1998 The image was generated using the second reprocessing algorithms and uses data that is excluded y 
the standard level -2 processing: atmospheric correction algorithm failure, land, sun glint, high total radiance, 
and clouds. In addition, stray light pixels are excluded. White areas indicate no data— continental land masses 
make up a large part of this region-while light gray indicates data present with occurrences of negative L wn 
values of less than 50%. The regions shaded black all have more than 50% occurrence of negative L WN values, 
indicating areas that are severely affected by negative Lwn- 


3.2.4. 3 Distribution of Negative Lvwv Values 

The occurrence of negative Lwn values, particularly 
in bands 1 and 2, has always been a problem in the Sea- 
WiFS processing scheme. Because the second reprocess- 
ing retained the actual values of the negative normalized 
water-leaving radiance instead of setting them to zero, it 
was possible to look at the distribution of the negative 
Lwn values. The geographic distribution of negative L W N 
values in the 412 nm band is presented in Fig. 6 for all 
level-2 data not masked by the level -2 flag conditions in 
an eight day period from 12-19 July 1998. Regions that 
had occurrences of negative Lwn values greater than 50% 
are shaded black. 

The negative Lwn distribution shows the largest prob- 
Jems to be in the coastal margins and in the extreme north- 
ern and southern latitude regions. The coastal waters have 
higher chlorophyll concentrations and are the most likely 
to have problems with turbid water or dust in the atmo- 
sphere affecting the aerosol correction. The non-coastal 


areas in the Southern Ocean (and to some degree in the 
extreme North Atlantic) with negative L W n values also 
have relatively high chlorophyll values, but should not have 
other coastal problems, such as continental aerosols. The 
amount of negative Lwn in these areas decreases as wave- 
length increases (not shown) and was fairly negligible at, 
and above, the 490 nm band. 

The distribution of negative Lwn values in higher lat- 
itude regions suggests the problem may be related to poor 
treatment of the atmospheric corrections at high solar ze- 
nith angles. The 12-19 July data were taken when the 
sun was at a latitude of 20°N. In that data (Fig. 6), the 
significant increase in negative Lwn values appeared to oc- 
cur at 40° S and possibly 80°N, which translates into solar 
zenith angles of around 60°. When the percentage of nega- 
tive Lwn values at 412 nm is plotted as a function of solar 
zenith angle (Sect. 3.3.6.2), a large increase in the percent- 
age of negative Lwn values is seen at a solar zenith angle 
of 60°. Similar behavior occurred for a period in January 
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1998 (not shown). This behavior may indicate a break- 
down in the plane parallel atmosphere assumption used to 
derive aerosol model behavior (Ding and Gordon 1994). 
This possibility was investigated further in preparation for 
the third reprocessing. 

An additional occurrence of negative L W n values was 
scattered throughout the open ocean areas, but unlike the 
coastal areas, the frequency of negative L WN values in 
these areas was relatively small. In the 510 and 555 nm 
bands (not shown), the distribution of open ocean negative 
Lwn values increased somewhat. This effect was primar- 
ily due to the influence of cloud shadows, noise, and stray 
light on the radiances. When the binning flags, which in- 
clude masks for low Lww(555) values, were included, this 
problem was almost entirely removed in band 1 and the 
other bands. 

The binning masks appeared to aid in removing some of 
the negative Lwn values, but large coastal areas remained 
severely affected by this problem. 

3.3 THE THIRD REPROCESSING 

The second reprocessing was able to improve on the al- 
gorithms used to perform the calibration and processing of 
the SeaWiFS data to water-leaving radiances and derived 
geophysical quantities. It also left many important ques- 
tions to be resolved. The most important concern was to 
determine what was causing the large number of negative 
Lwn values in coastal areas. The effort to understand and 
solve this problem and the algorithms chosen for the third 
reprocessing are presented here. 

In the period leading up to the third reprocessing, a 
number of improvements were also made in the calibra- 
tion and the atmospheric correction algorithms. Calibra- 
tion improvements were possible with more measurements, 
which led to a greater understanding of instrument perfor- 
mance. Another significant change was the migration of 
the level -2 processing algorithms to a new program, MSI 12, 
which greatly improved the ease of use and the incorpora- 
tion of new algorithms. 

The MSI 12 program permitted the selection of many 
more intermediate and final parameters. With this mi- 
gration, several basic changes were made in the default 
operational products and in the available flags. The major 
operational changes are summarized in Table 4 and brief 
descriptions of the changes are presented in the following 
sections along with the results of the third reprocessing 
(Sect. 3.3.6). 

3.3.1 Calibration Related Changes 

The changes to the calibration included the routine 
characterization of the detector degradation, as well as, 
a better understanding of instrument behavior. Increased 
use of the SeaWiFS data for cloud and aerosol studies re- 
vealed other calibration problems, which were solved for 
the third reprocessing. 


3.3. 1.1 Expanded Time- Dependent Calibration 

The SeaWiFS CVT has updated the time-dependent 
calibration as new lunar calibration measurements become 
available each month. Previously (Sect. 3.2. 1.2), a time de- 
pendence could only be seen in the two near-infrared bands 
(765 and 865 nm). A longer time period of observations, 
and a more accurate treatment of the radiance model of 
the moon, made it possible to detect and implement a cor- 
rection for time-dependent changes in the radiometric re- 
sponse in many of the remaining bands (Eplee and Barnes 
2000 ). 

The largest changes now seen in bands 7 and 8 are a 
2% and 10% decrease, respectively, while bands 1 and 6 
(412 and 670 nm) show 1% decreases, and bands 2 and 5 
(443 and 555 nm) show 0.5% decreases. Bands 3 and 4 (490 
and 510 nm) are assumed to have no change. The SeaWiFS 
calibration now contains time-dependent corrections for all 
these bands. 

3.3. 1.2 Bilinear Gain Adjustments 

During work to derive a measure of absorbing aerosols 
(Hsu et al. 2000 and Fukushima et al. 1999), it was de- 
termined that there was a discontinuity in the frequency 
distribution of total radiances which occurred at the region 
where the instrument gain changed to a lower sensitivity 
(the so-called knee point). This indicated that the relation- 
ship between raw satellite counts and total radiance was 
incorrect for radiance values above the knee. The radiance 
value at which this discontinuity occurred was above the 
range where most ocean color processing is performed and, 
therefore, has no effect on the oceanic products. It does, 
however, have an effect on the analysis of dense aerosols, 
absorbing aerosols, and clouds. It was assumed that the 
laboratory setting of the knee point was either incorrect or 
had shifted. 

The knee points were redetermined by finding the set- 
tings that minimized the discontinuity in the frequency 
distributions of the total radiances (Eplee and Patt 2000). 
The implementation of this change reduced the radiances 
above the knee by 0.8% on average, and increased the ra- 
diance below the knee by about 0.1%. 

3.3. 1.3 Revised Vicarious Calibration 

The calibration of bands 1-7 was performed for the 
third reprocessing in much the same way as it was done for 
the second reprocessing (Sect. 3.2. 1.3, Eplee and McClain 
2000a, and Robinson and Wang 2000). This calibration 
is performed every time the level -2 processing algorithms 
are modified to ensure that the water-leaving radiances 
produced by the algorithms match the ground-truth mea- 
surements from MOBY. 

The vicarious calibration is also performed to include 
additional SeaWiFS and MOBY measurements taken over 
time. 
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Table 4. Major changes and their observable effects in the third reprocessing. 


Change 

Observable Effect 

Calibration: 


Expanded time-dependent calibration 
Bilinear gain adjustments 
Revised vicarious calibration 

Better temporal stability 

Smoothed transition to higher radiances 

None 

Temperature correction 

None 

Atmospheric Correction: 

Modified aerosol model selection 
Epsilon value extrapolation 
Improved Rayleigh computation 
Pressure dependent transmittance 
Near-infrared Lwn adjustment 

Reduced discontinuity in some fields 
Better retrieval of low chlorophyll values 

Significant decrease in negative Lwn a t high solar zenith angles 
Small Lwn changes 

Significant decrease in negative Lwn an d high chlorophyll 
values in coastal areas 

Other: 

Modified C a 'K computation 
Spectral-dependent whitecap correction 
Sun glint correction 

Eliminated bad data values in binned data 
Small reductions in negative Lwn values 
Greatly improved aerosol optical thickness 

Ozone data interpolation 
Improved A(490) algorithm 
Absorbing aerosol flagging 
Trichodesmium flagging 
Improved chlorophyll algorithm 
Out-of-band correction 
Navigation Improvements 
Improved coastal data inclusion 

None 

Better K( 490) values in turbid water 
Reduced binning of contaminated data 
None 

Slight lowering of low chlorophyll values 
Better compatibility to in situ measurements 
Improved location of data 
More coastal retrievals 

Structural: 

Product and flag updates 
Level -2 program code changes 

Better flag specificity 

Greater flexibility and faster updates 


3 .3. 1.4 Temperature Correction 

A correction was made to the level -0 to -1 conversion 
software to correctly unpack the focal plane temperatures 
used in the calibration process. The effect of this correc- 
tion on the water- leaving radiance values was less than 
0.5% in bands 1-4, 1.5% in band 5, and 3.0% in band 6. 
Recalibration will remove most of these changes. 

3.3.2 Atmospheric Correction Changes 

A number of improvements were made in the algo- 
rithms that remove the radiance contributions of the at- 
mosphere. As with the changes in the calibration (Sect. 
3.3.1), some of the changes were made to correct artifacts 
seen in the data, and others were made to perform a more 
accurate atmospheric correction. 

3.3.2. 1 Modified Aerosol Model Selection 

It was discovered that the aerosol radiance fields con- 
tained noticeable discontinuities in open ocean areas. The 
discontinuities, which paralleled lines of constant scatter- 
ing angle, propagated noticeable and unwanted artifacts 


in the SeaWiFS products. The effect became more notice- 
able after smoothing of the near-infrared radiance fields a 
method under investigation for improving the data quality 
(e.g., to reduce speckling). 

The discontinuity was caused by the transition between 
the oceanic aerosol model with 90% humidity and other 
models. This model was removed, leaving the oceanic 
aerosol model with 99% humidity to handle the very clear, 
oceanic aerosol conditions. To preserve the 12 model set, 
the coastal model with 70% humidity was reinstated to the 
model suite. 

3. 3. 2. 2 Epsilon Value Extrapolation 

Even with the oceanic aerosol models that extend the 
coverage to a larger aerosol type range, there were many 
occasions when this range was exceeded. On these occa- 
sions, the e value, which is an indicator of the aerosol type, 
has a value lower than that of the lowest aerosol model. In 
the second reprocessing, the e value was used to interpolate 
between two aerosol models, or in the case of an e value 
below the lowest model, the value of the lowest model was 
used (instead of the true value). The result was the re- 
moval of more aerosol radiance than if the e value was 
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extrapolated. The change made for the third reprocessing 
was to extrapolate the e values for the other bands using 
an analytical function, so a better aerosol correction could 
be made. 

For test cases in open ocean areas where e value ex- 
trapolation would occur, up to 24% of the Lwn retrievals 
benefitted from the extrapolation. No significant changes 
were observed in the chlorophyll a fields, but the 412 nm 
water-leaving radiances were increased by an average of 

0.1 mW cm” 2 /im ” 1 sr” 1 , an increase of about 6%. 

3. 3. 2. 3 Improved Rayleigh Computation 

The Rayleigh radiance algorithm, which was used for 
the second reprocessing, contained no correction for ocean 
surface roughness. The effect of ocean roughness is neg- 
ligible at low and moderate solar zenith angles, but be- 
comes significant at solar zenith angles greater than 60°. 
Without an adjustment for wind roughening, the Rayleigh 
radiance estimate is too high at large solar zenith angles 
and greater than zero wind speeds. This results in a reduc- 
tion of the Lwn values and an increase in the occurrence 
of negative Lwn values. Tests of level -2 retrievals with 
the wind-dependent Rayleigh algorithm showed that neg- 
ative water-leaving radiances were reduced in the 412 nm 
band from 80% to 40% at solar zenith angles of 70°, where 
the problem is most severe. The Lwn values at high solar 
zenith angles also increased significantly. The effects are 
addressed further in Sect. 3. 3. 6. 2. 

3. 3. 2. 4 Pressure Dependent Transmittance 

Prior to the third reprocessing, the diffuse transmit- 
tance used in the computation of the normalized water- 
leaving radiance was derived from a look-up table with 
dependence only on the aerosol model and the viewing 
geometry. Currently, a pressure correction (through the 
Rayleigh optical thickness) is also applied to both the sun- 
to-surface and surface-to-satellite diffuse transmittance 
(Wang 1999b). The effect on Lwn values is small: a 
large change in the atmospheric pressure of 30 mbar causes 
about a 1% change in Lwn values at 412 nm, while at 
865 nm, the change is only 0.05%. 

3.3.3 Near-Infrared Lvvtv Value Adjustment 

One of the problems noted in the second reprocessing, 
which was examined closely in preparation for the third 
reprocessing, was the occurrence of negative water-leaving 
radiances in the SeaWiFS bands. It was apparent that 
especially in coastal regions, the water-leaving radiances in 
the 412 and 443 nm bands were negative for two reasons: 

a) High chlorophyll content, which depressed the blue 
radiances, made water-leaving radiances approxi- 
mately zero; and 

b) Turbid coastal waters confuse the aerosol correc- 
tion (absorbing aerosols could also be responsible), 


resulting in the overestimation of aerosol radiance, 
especially in the blue bands, thereby decreasing the 
already low radiance value, frequently causing it to 
become negative. 

The negative water- leaving radiances were obviously unre- 
alistic, although radiances approaching zero were common 
in turbid waters for bands 1 and 2. The negative Lwn 
values made the affected pixels useless for algorithms that 
rely on the blue bands as input. In the same regions having 
the negative Lwn problem, anomolously high chlorophyll 
was found, which again indicated that the blue bands were 
being excessively depressed. 

The ocean science community responded to this prob- 
lem with a number of algorithms to explain the problem 
and to get better Lwn estimates. The solutions fell into 
two categories: 

1. The assumption of zero water-leaving radiance in 
the near-infrared bands at 765 and 865 nm is wrong 
in high chlorophyll or turbid waters. The algo- 
rithms attempt to use the chlorophyll values to es- 
timate the water-leaving radiances in these bands; 
the algorithms were contributed by D. Siegel (Siegel 
et al. 2000), R. Arnone, and R. Stumpf. 

2. After the near-infrared contribution is taken into 
consideration, an assumption can be made about 
the value of the 412 nm water-leaving radiance so 
that it and the other Lw values do not become neg- 
ative (R. Stumpf, pers. comm.). 

The most promising combinations of these methods were 
tested, including the use of band 6 (instead of band 7) with 
band 8 in the aerosol determination. 

The methods were examined in detail for a number 
of test scenes that were affected by the problem. In ad- 
dition, the methods were tested in general on two eight- 
day periods of SeaWiFS GAC data — one in January 1998 
and one in July 1998. The results were then compared to 
the method used in the second reprocessing. Although in- 
vestigation of the other methods will continue, the Siegel 
method (Siegel et al. 2000) was chosen for the third repro- 
cessing, because of its simplicity, and because it produced 
significant decreases in the amount of negative Lwn values 
while lowering coastal chlorophyll values to more reason- 
able levels. 

3.3.4 Other Changes 

Several other changes were made in the level -2 and -3 
processing in preparation for the third reprocessing. As a 
result of the level -2 data set format change (Sect. 3.3.5), 
two new flagging algorithms were added to indicate the 
presence of absorbing aerosols and Trichodesmium . Im- 
provements were made to the chlorophyll, diffuse atten- 
uation, ozone, and whitecap algorithms. A glint correc- 
tion scheme yielded better atmospheric optical depths and 
an out-of-band correction permitted better comparison of 
Lwn with in situ values. 
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3.3.4. 1 Modified C a :K Computation 

After the second reprocessing was completed, it was 
discovered that some level -3 binned products had infinite 
values stored in some bins for the parameter containing the 
ratio of the chlorophyll a concentration values to the diffuse 
attenuation coefficient at 490 nm, C a :A(490), hereafter re- 
ferred to as C a -K. The problem was introduced because 
of the new binning policy, which allowed the binning of 
negative L W n data, combined with some shortcomings in 
the binning code. 

The algorithm for calculating A (490) in the level -2 pro- 
cessing was set to return a zero value for A (490) if the 
Lwn at either 490 or 555 nm was less than zero. When a 
A(490) value of zero was used in the binner during the cal- 
culation of C a -K, a value of infinity was assigned to that 
bin. This problem was compounded for products covering 
longer time periods because once a bin had an infinite value 
assigned to it, no amount of averaging would lessen or re- 
move it; therefore, eight-day, monthly, and yearly binned 
products would accumulate more bins with infinite C a A 
values. This problem was fixed in both the A(490) algo- 
rithm in the level -2 processing code and in the C a K com- 
putation in the binner. In the A (490) algorithm, the com- 
putation was modified so that zero A (490) values would 
not be produced. 

The following rules were designed to deal with the Lwn 
conditions: 

a. If L wn (490) and Lwn (555) > 0- then compute 
A(490) normally. 

b. If Lwn (490) < 0, then A(490) = 6.4 m- 1 . 

c. If Lwn ( 490) > 0 and L wn at 555 nm < 0, then 
A (490) = 0.016 m- 1 . 

The actions taken for zero Lwn values represent the max- 
imum A(490) value that can be stored in the level-2 data 
set and the lowest possible A(490) value which the algo- 
rithm can produce. Tests with SeaWiFS data showed that 
the default values agreed well with the surrounding A (490) 
values. 

The algorithm for computing C a -K in the binner was 
also adjusted to use a A(490) value of 0.016 m- 1 in the 
calculation if a value of zero was encountered. In addition, 
a modification was made in the masking conditions in the 
event of an atmospheric failure so that fewer high A (490) 
values, erroneously derived over open ocean areas, would 
be binned. If the L w value for any of the 490, 510, or 
555 nm bands are less than zero, the atmospheric warning 
flag is set in the level -2 file and applied as a mask in the 
level -3 binning. 

3.3.4. 2 Changes to Whitecap Correction 

For the third reprocessing, three changes were made 
in the whitecap correction. The first change was to use a 
correction with a spectral dependence (Frouin et al. 1996). 


The second change was to increase the strength of the cor- 
rection by approximately 65% from the values used in the 
second reprocessing. The third change was to limit the 
whitecap correction for wind speeds above 8ms to the 
value found at 8ms- 1 . This limit was set to avoid over- 
corrections for whitecaps. This strategy incorporated more 
recent information on whitecaps (Moore et al. 2000), but 
avoided overcorrection at high wind speeds. 

The effects on the radiances were small, but noticeable. 
No significant changes were observed for areas with wind 
speeds less than 10 m s' 1 . For areas with wind speeds from 
10-15 ms -1 , the occurrence of negative water-leaving radi- 
ances at 412 nm was reduced by 3%; and for wind speeds 
above 15ms -1 , negative Lwn values at 412 nm were re- 
duced by 10%. 

3.3. 4.3 Sun Glint Correction 

The SeaWiFS Project found that the aerosol optical 
thickness was noticably higher for areas surrounding the 
subsolar point. The atmospheric correction was account- 
ing for glint outside the glint mask as additional aerosol 
radiance. This did not noticably affect the Lwn and chlo- 
rophyll a retrievals, but it inflated the aerosol optical thick- 
ness, thereby making it less useful. 

The level-2 processing program already calculated an 
estimate of the glint radiance, so the program was changed 
to remove the glint radiance outside the glint mask as a 
part of the processing. The glint removal was made more 
robust by performing a second iteration with a better value 
of the aerosol optical thickness. Details of the glint correc- 
tion are presented in Wang and Bailey (2000). 

3,3. 4 . 4 Ozone Data Interpolation 

In previous reprocessings, TOMS ozone values (the pri- 
mary SeaWiFS ozone source) were interpolated to the time 
of the SeaWiFS data pass assuming that the TOMS data 
were all taken at 12:00 UTC. In fact, the TOMS data 
were collected over the daylight side of the Earth from 
east to west with a phasing very close to SeaWiFS (the 
data were taken almost at the same time). The ancillary 
data selection and interpolation routines in the level -2 pro- 
cessing were modified to use this information to calculate 
ozone fields, which would be more representative of the 
actual conditions at the time of the pass (Ainsworth and 
Patt 2000). The effect is the smallest around 12:00 UTC, 
which is near the prime meridian, and the greatest around 
0:00 UTC near ±180° longitude. 

3. 3. 4- 5 Improved K(490) Algorithm 

An improved (490) algorithm (Mueller 2000) was im- 
plemented for the third reprocessing. The previous algo- 
rithm used Lw (443) and had errors in regions of highly 
turbid water and in bloom situations where L w (443) can 
be underestimated. The new algorithm estimates K(490) 
by using the ratio of Lwjv(490):Lwn(555). 
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3. 3. 4-6 Absorbing Aerosol Flag 

In many oceanic regions, and especially off the west 
coast of Africa, absorbing aerosols (suspended dust) are 
not handled by the aerosol determination algorithm. Cur- 
rently, the suite of aerosol models do not contain absorbing 
aerosol models. For the third reprocessing, an algorithm 
was implemented to detect significant amounts of absorb- 
ing aerosols (Hsu et al. 2000). A new level-2 flag was 
made to indicate measurements in regions of excessive ab- 
sorbing aerosols. This flag is used to exclude data during 
the level -3 binning phase. 

3.3,4* 7 TVichodesmium Flag 

An algorithm was implemented in the third reprocess- 
ing that detects the presence of Trtchodesmium bloom con- 
ditions (Subramaniam et al. 2000). The existence of a 
bloom is indicated by one of the new level -2 flags, which 
currently is not used for any level ~3 exclusion. 

3.3.4- & Improved Chlorophyll a Algorithm 

The chlorophyll a algorithm used for the third repro- 
cessing, the so-called ocean chlorophyll 4 (OC4) algorithm 
(O'Reilly et al. 1998), was used in place of the ocean chlo- 
rophyll 2 (OC2) algorithm (second reprocessing), A larger 
in situ data set was used to refine the OC4 algorithm 
(O’Reilly et al. 2000), and it also has the property of 
not reporting negative chlorophyll a values as the previ- 
ous OC2 algorithm did. The range of chlorophyll in the 
output is expanded at the low end to include values from 
0-0.009 mg m -3 in steps of 0.001. Previously, the lowest 
reported chlorophyll value was 0.01 mg m~ 3 . In general, 
the new algorithm was found to slightly reduce the chlo- 
rophyll a values in the range from 0.01-0.05 mg m“ 3 . 

3. 3.4- 9 Out-of-Band Correction 

The eight SeaWiFS bands have broad response func- 
tions compared to most instruments that measure L w in 
the field. Although the processing for SeaWiFS water- 
leaving radiances includes the out-of-band response, the 
resulting Lw values retain the initial broad response func- 
tion. For the third reprocessing, an out-of-band correction 
was applied to the Lw values as the default (Wang et al. 
2000). The correction has the greatest effect at 555 nm 
in low chlorophyll concentrations (high blue radiances), 
where the corrected L W n values can be 5-10% lower than 
the uncorrected values. The use of the out-of-band cor- 
rection slightly decreased the chlorophyll values by about 
15% in the range from 0.01-0.05 mgm -3 . 

3.3.4- 10 Navigation Improvements 

Improvements were made to the navigation algorithms 
in the level -0 to -1A software to reduce the seasonal vari- 
ations in geolocation accuracy and to handle operational 
changes in the available navigation data from the satellite. 


3.3-4-11 Improved Coastal Data Inclusion 

Some instances were observed when the radiance at 
412 nm— after being corrected for Rayleigh, glint, ozone, 
and whitecaps— was negative for otherwise clear regions. 
This condition occurred infrequently along coastal areas, 
but affected significant regions in any one single data pass. 
The atmospheric failure masking used in the second re- 
processing used this test to exclude observations that had 
these conditions in any of the eight SeaWiFS bands. In 
the third reprocessing, this exclusion was removed for the 
412 nm band so more coastal areas could be processed. 
This change was observed to increase the number of re- 
trievels at many coastal regions for approximately one out 
of every four passes. 

3.3.5 Major Product and Flag Changes 

Since the SeaWiFS launch, the Project has been cre- 
ating the same suite of products and flags (McClain et al. 
1995 and McClain 2000). With the third reprocessing, the 
Project took the opportunity to revise the product suite 
and to expand the flag set to respond to new demands. 

3.3.5. 1 Product Suite Changes 

Table 5 shows the operational product suite used in the 
third reprocessing and the product names used in SeaDAS. 
The Coastal Zone Color Scanner (CZCS) pigment product 
(CZCS_pigment) was removed because there was little de- 
mand for it, and because a simple equation can be used to 
derive the pigment from the chlorophyll a product. In ad- 
dition, the two products of the aerosol radiance at 670 nm 
and at 865 nm (La_670 and La_865, respectively) were re- 
moved in favor of the normalized water-leaving radiance 
at 670 nm (nLw_670). The Angstrom coefficient at 510 nm 
(angstrom_510), was included because many researchers 
use the Angstrom coefficient to characterize the aerosol 
type. 

Table 5. Operational products for the third repro- 


cessing. 

Product Name 

Description 

nLw_412 

Lwn values at 412 nm 

nLv-443 

Lwn values at 443 nm 

nLw_490 

Lwn values at 490 nm 

nLw.510 

Lwn values at 510 nm 

nLw_555 

Lwn values at 555 nm 

nLw-670 

Lwn values at 670 nm 

angstrom_510 

Angstrom coefficient at 510 


and 865 nm 

chlor_a 

Chlorophyll a concentration 

K_490 

Diffuse attenuation coefficient 


at 490 nm 

eps_78 

Epsilon at 765 and 865 nm 

tau_865 

Aerosol optical thickness at 


in 865 nm 
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Table 6. Flags for the third reprocessing. The “Flag Name” column denotes the flag names as of the third 
renrocessine whereas the “Old Name” column is the flag name used m the second reprocessing. The Mask 
In” columns’ indicate that no geophysical data is created in the level -2 (L2) or level -3 (L3) data set if the flag 


conditions 

Flag 

Number 

marKea i ex 

Flag 

Name 

1 st ror wi»l 
Mask In 
L2 L3 

UUSCl VOLIKJli. 

Old 

Name 

Description 

1 

ATMFAIL 

Y 

Y 

EPS I LON 1 

Atmospheric algorithm failure 

2 

LAND 

Y 

Y 

LAND1 

Land 

3 

BADANC 



ANCIL1 

Missing ancillary data 

4 

HIGLINT 

Y 

Y 

SUNGLINT1 

Sun glint contamination 

5 

HILT 

Y 

Y 

HIGHLT1 

Total radiance above the knee in any band 

6 

HISATZEN 


Y 

SATZEN1 

Satellite zenith angle above the limit 

7 

C0ASTZ 



C0ASTZ1 

Shallow water 

8 

NEGLW 



NEGLW 1 

Negative water-leaving radiance in any band 

9 

STRAYLIGHT 


Y 

STRAYLIGHTl 

Stray light contamination 

10 

CLDICE 

Y 

Y 

CLDICE1 

Clouds or ice 

11 

C0CC0LITH 


Y 

C0CC0LITH1 

Coccolithophore bloom 

12 

TURBIDW 


Y 

TURBIDV1 

Turbid (Case-2) water 

13 

HIS0LZEN 


Y 

S0LZEN1 

Solar zenith angle above the limit 

14 

HITAU 



HIGHTAU1 

High aerosol concentration 

15 

L0WLW 


Y 

L0WLW1 

Low water-leaving radiance at 555 nm 

16 

CHLFAIL 

t 

Y 

CHL0R1 

Chlorophyll not calculable 

17 

NAVWARN 


Y 


Questionable navigation (tilt change) 

18 

ABSAER 


Y 


Absorbing aerosol index above the threshold 

IQ 

TRICH0 




Trichodesmium bloom condition 

j. j 
20 

X lVlvllv 

MAXAERITER 


Y 


Maximum number of iterations in the NIR algorithm 

21 

M0DGLINT 




Glint corrected measurement 

22 

CHLWARN 




Chlorophyll is out of range 

23 

ATMWARN 


Y 


Epsilon value is outside the reasonable range or Lw 






in 510, at 490, or 555 nm is less than zero 

24 

DARKPIXEL 




Rayleigh corrected radiance is less than zero for 






any band 

_ A 


■(■ The chlorophyll value is not computed, but first 


Similar changes were made in the level _ 3 binned prod- 
ucts. In place of the La_670 and CZCS.pigment products 
are the nLw_670 and angstrom_510 products. 

3. 3. 5. 2 Flag Changes 

For previous reprocessings, the suite of flags was limited 
to 16 by the format of the level -2 data set. In the third 
reprocessing, the available room for flags was expanded to 
32, of which only 24 flags are currently defined. Table 6 
lists these flags and their status as masks in excluding data 
in the level-2 and - 3 operational products. 

The original 16 flags have much the same meaning as 
before with some exceptions. In the second reprocess- 
ing, the atmospheric algorithm failure flag (EPSIL0N1) in- 
cluded many conditions that prevented the calculation of 
the aerosol radiances and, thus, the Lwn values. It also 
indicated times when the e value was outside reasonable 
limits defined by the standard aerosol models; hence, the 
name EPSILON 1. Now, the atmospheric algorithm failure 
flag (renamed to ATMFAIL), only indicates conditions where 


Lwn values could not be calculated. The new ATMWARN flag 
indicates when the t value is outside the reasonable range, 
but Lwn values could be calculated. It also indicates ob- 
servations where L w values in the 490, 510, and 555 nm 
bands are negative. A serious problem condition, when 
the Rayleigh corrected radiances are negative, is indicated 
in the new DARKPIXEL flag so that this condition can be 
monitored more easily. 

In the second reprocessing, the chlorophyll algorithm 
failure flag (CHL0R1) indicated two problems: a) chloro- 
phyll could not be computed, and b) chlorophyll exceeded 
the high threshold. Now, the CHLFAIL flag only indicates 
when chlorophyll can not be calculated, because the in- 
put Lwn values are less than zero, or when the calcu- 
lated chlorophyll is outside physical limits (greater than 
640 mg nT 3 ). The new CHLWARN flag signals chlorophyll 
values which exceed the high value that can be stored in 
the product (greater than 64mgm' 3 ), or very small values 
(less than 0.01 mgm -3 ). 

The remaining six new flags address new conditions or 
are used to more consistently handle current conditions. 
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Table/. Average chlorophyll in the Chesapeake Bay derived from in situ measurements (SeaWiFS second and 
third data reprocessings). The mean chlorophyll a values (C a ) are in milligrams per cubic meter (mgm -3 ) The 
standard deviation ( cr ) and number of observations are also given. 


Date 

In Situ 

C a cr Obs. 

Second Reprocessing 
C a cr Obs. 

Third Reprocessing 
C a a Obs. 

11-19 April 1998 
4-12 August 1998 
19-23 October 1998 

11.93 11.31 91 

10.83 10.61 89 

7.43 4.08 67 

42.64 18.42 1755 

25.93 13.93 1804 

16.85 9.76 1815 

16.01 7.96 1960 

11.06 3.71 1810 

7.96 2.97 2133 


The NAVWARN flag is used primarily to indicate where less versus the second reprocessing. For a number of chloro- 
reliable navigation is expected, such as when the instru- phyll ranges in the second reprocessing, the statistics are 
ment tilt is changing. The level-3 binning can use the plotted with the matching third reprocessing chlorophyll 
NAVWARN flag to decide on binning instead of having to use values. In the 0. 2-1.0 mgm -3 range, the third reprocess- 
the tilt indicator in the level-2 data. The ABSAER and ing shows a slight increase in the chlorophyll value of about 
TRICHO flags indicate the existence of absorbing aerosols 10%, as was seen in the histograms (Fig. 7). For chloro- 
and Trichodesmium blooms detected by new algorithms, phyfl values greater than 2 mgm -3 , the third reprocess- 
respectively. The MAXAERITER flag indicates when the near- ' n § significantly lowered the chlorophyll value. In the 5- 


infrared (NIR) correction algorithm (Sect. 3.3.3) has been 
unable to converge on a NIR L WN estimate. Finally, the 
MODGLINT flag indicates pixels where the glint correction 
was applied. 

3.3.6 Third Reprocessing Analysis 

The software for the third reprocessing was updated 
with all the changes mentioned in Sects. 3.3. 1-3. 3. 5. The 
vicarious calibration for bands 1-7 was performed with the 
new algorithms. The new algorithms and the new calibra- 
tion were used to process two eight-day periods: 17-24 
January 1998 and 12-19 July 1998. The results of these 
runs, in comparison to the second reprocessing runs, is 
discussed next. 

3.3.6. 1 Chlorophyll a 

Figure 7 is a plot of the global distribution of chlo- 
rophyll a concentration for the January and July periods. 
The inclusion of binned chlorophyll in the 0-0.009 mgm -3 
range for the third reprocessing can be seen, as can the 
artificial accumulation of values, at 0.01 mgm -3 , for the 
second reprocessing. The third reprocessing produced a 
greater number of chlorophyll values in the 0-0.04 mgm -3 
range than the second reprocessing, which indicates that 
the chlorophyll values in this range are reduced. For a test 
area in the Pacific located in one of the lowest chlorophyll 
areas in the data, around (135°E,12°N), the average chlo- 
rophyll value dropped from 0.057-0.0381 mgm -3 , a 34% 
decrease. In the chlorophyll range above 0.07 mgm -3 , the 
third reprocessing increased the chlorophyll values. The 
magnitude of the changes in July period are less, but still 
follow the same trend. 

Figure 8 shows the change that occurred in chloro- 
phyll a concentration more clearly for the intermediate and 
high values of chlorophyll, and it summarizes a scatterplot 
comparing binned chlorophyll from the third reprocessing 


10 mgm -3 range, the third reprocessing lowered the chlo- 
rophyll values by about 40%, and lowers it even more for 
higher values. This lowering of high chlorophyll is primar- 
ily due to the Siegel NIR algorithm (Siegel et al. 2000) 
and greatly reduces chlorophyll values in all the coastal 
and upwelling regions. 

To illustrate how the third reprocessing improves chlo- 
rophyll a retrievals in coastal waters, SeaWiFS chlorophyll 
values were compared with in situ pigment data from the 
Chesapeake Bay. Three time periods in 1998 were selected 
for this analysis: 11-19 April, 4-12 August, and 19-23 Oc- 
tober. For each time period, HRPT data were processed 
to level -2 chlorophyll using both the second and third re- 
processings and methodologies, and were then space- and 
time-binned. In situ data (originally provided by L. Hard- 
ing, University of Maryland), were obtained from SeaBASS 
(see Sect. 7.1.1 for an overview of SeaBASS). Only the 
lower Chesapeake Bay, south of 38.5°N latitude, was con- 
sidered in the analysis. 

For all three time periods, the mean SeaWiFS-retrieved 
chlorophyll value exceeds that of the mean in situ value. 
The third reprocessing values, however, are 62.5, 57.3, and 
52.8% lower than the second reprocessing values, for April, 
August, and October, respectively (Table 7). Additionally! 
the ratio of mean SeaWiFS-to-mean in situ chlorophyll val- 
ues is reduced from 3.57-1.34 for April, 2.39-1.02 for Au- 
gust, and 2.27-1.07 for October. Results from these com- 
parisons (especially the latter) indicate that the method- 
ology for the third reprocessing significantly improves Sea- 
WiFS chlorophyll a retrievals in coastal regions. The num- 
ber of successful chlorophyll retrievals also increased by 
more than 10% for two of the three time periods. 

3. 3. 6. 2 Negative Values 

The use of the Siegel NIR algorithm (Siegel et al. 2000) 
and the wind speed correction to the Rayleigh radiance al- 
gorithm (Sect. 3. 3. 2.3) significantly helped reduce the num- 
ber of negative L WN values globally. In the January time 
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Fig. 7. The SeaWiFS level -3 global chlorophyll a frequency distribution for the eight day period of: a) 
17-24 January 1998, and b) 12-19 July 1998, made with the second reprocessing (dotted curve) and the third 
reprocessing (solid curve). 



Fig. 8. Comparison of the chlorophyll retrieved in the third reprocessing (y-axis) to that retrieved in the 
second reprocessing (x-axis) for: a) 17-24 January 1998, and b) 12-19 July 1998. For level -3 binned chlo- 
rophyll a concentration from the second reprocessing in 11 ranges from 0.02-64.0 mg m , the median third 
reprocessing chlorophyll value is plotted as diamonds. The nearest set of bars is the range of data in the 
second-to-third quartiles, while the outer bars show the extreme range of the third reprocessing chlorophyll 
values. Chlorophyll values greater than 2mgm -3 were significantly reduced in the third reprocessing. 
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Fig. 9. The percent occurrence of negative L W n values at 412 nm for the eight-day period from 12-19 
July 1998. This image was generated using the third reprocessing algorithms and used data excluded by the 
standard level-2 processing: atmospheric correction algorithm failure, land, sun glint, high total radiance, and 
clouds. In addition, stray light pixels were excluded. White areas indicate no data— continental land masses 
make up a great part of this region— while light gray indicates data present with occurrences of negative 
L W n values of less than 50%. The black regions all have more than 50% occurrence of negative Lwn values, 
indicating areas that were severely affected by negative Lwn values. In regions where the solar zenith angle 
is high, such as the Southern Ocean and along the coast, the amount of negative Lwn values were reduced 
relative to the second reprocessing (Fig. 6). 


period, the number of negative L W n values in the binned 
data for the 412 nm band was reduced from 2.5-2. 1% of 
the total, while in the July time period, the reduction was 
from 5. 2-4. 7%. Figure 9 shows the geographical regions, 
in the July time period, wherein the percentage of negative 
Lwn values are higher than 50% for the third reprocess- 
ing. Compared with the same display for the second re- 
processing (Fig. 6), the reduction in negative Lwn values 
at high solar zenith angles, and especially in the South- 
ern Ocean, is dramatic. The improved performance of the 
atmospheric correction with the wind-dependent Rayleigh 
radiances is the main reason why the solar zenith limit for 
useful data was increased from 70-75°. Slight reductions 
can also be seen in the occurrence of negative Lwn values 
along coastal areas. 


When the wind speed correction to the Rayleigh radi- 
ances is made, the amount of negative Lwn values is re- 
duced. This effect is illustrated in Fig. 10, which shows the 
percentage of negative Lwn as a function of solar zenith 
angle for the 412, 490, and 555 nm SeaWiFS bands in the 
July time period (the other bands and the January time 
period show similar behavior). In the 412 nm band, at 
solar zenith angles from 65-70°, the percentage of mea- 
surements with negative Lwn values dropped from 60% to 
20%, with similar decreases in the other bands. 

Figure 11 shows the mean of the Lwn values as a func- 
tion of solar zenith angle for the same SeaWiFS bands 
in the July time period. The decrease in the mean Lwn 
values, which begins at solar zenith angles of about 60°, 
could lower the Lwn values in the 412 nm band to less than 
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Fig. 10. The percentage of negative Lwn values as 
a function of the solar zenith angle for the eight-day 
period of 12-19 July 1998 for the second (dotted curves) 
and third (solid curves) reprocessings. The SeaWiFS 
bands at a) 412, b) 490, and c) 555 nm are shown. 



Fig. 11. Mean Lwn values as a function of solar zenith 
angle for the eight-day period of 12-19 July 1998 for the 
second (dotted curves) and third (solid curves) repro- 
cessings. The SeaWiFS bands at a) 412, b) 490, and 
c) 555 nm are shown. 
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zero in the second reprocessing (Fig. 11a). In addition, the 
555 nm band average (Fig. 11c) dropped below 0.15, which 
is a masking threshold for removing data contaminated 
by cloud shadows. The third reprocessing significantly in- 
creased these means at high solar zenith angles and greatly 
reduced the aforementioned problems in both the 412 and 
555 nm bands, as well as yielding increased Lwn values 
in the other bands and restoring a more constant value of 
Lwn in the green bands (i.e., 510 and 555 nm). 

In general, the number of level -3 bins filled in the eight- 
day time bins of the third reprocessing increased over those 
in the second reprocessing by 1.2% and 2.4% for the Jan- 
uary and July test periods, respectively. This corresponds 
to an increase of 1.8xl0 6 and 3.8xl0 6 km 2 in ocean cov- 
erage for the same respective time periods. Much of this 
increase is due to the improved, wind-dependent Rayleigh 
radiance algorithm, and the opening of the solar zenith 
angle cutoff for binning from 70-75°. 

More results of the quality of the SeaWiFS Lwn val- 
ues and derived chlorophyll a concentrations can be found 
in Bailey et al. (2000). These comparisons show marked 
improvement in the Lwn and chlorophyll values compared 
to in situ measurements. 

3*4 CONCLUSIONS 

The second reprocessing repaired many large problems 
and brought about favorable changes in both the chlo- 
rophyll and Lwn distributions. The changes increased 
very low chlorophyll values in open-ocean regions, and de- 
creased the values of high chlorophyll values greater than 
l.Omgm -3 and in many cases, allowed more chlorophyll 
retrievals to be made and binned in the level -3 data. Lwn 
values generally increased, although the prevalence of 


negative Lwn was recognized as a major problem remain- 
ing. Negative Lwn values were found to be most prevalent 
in bands 1 and 2, at higher solar zenith angles, and around 
coastal areas where turbid, high chlorophyll concentrations 
exist. 

The third reprocessing increased the usefulness of Sea- 
WiFS data in both oceanic and atmospheric applications. 
The negative Lwn problem, which was recognized and 
characterized in the second reprocessing, was addressed 
and significantly reduced in the third reprocessing through 
the use of a) a correction of the NIR L W n values, and b) 
a wind-dependent correction to the Rayleigh radiance cal- 
culation. Abnormally high coastal chlorophyll a concen- 
trations were decreased, and data at higher solar zenith 
angles were retrieved, as a result of these improvements. 
A chlorophyll algorithm that benefits from an increased 
observation set improved the retrieval of chlorophyll from 
the Lwn data. These and many other changes combined 
to improve the agreement between SeaWiFS L W n, chloro- 
phyll, and aerosol optical thickness with corresponding in 
situ observations. 

In summary, the changes made for the second and third 
reprocessings have been wide ranging and have had a pos- 
itive influence on all of the SeaWiFS products. Even so, 
there are many more possible improvements in future re- 
processings that will be explored. The occurrence of nega- 
tive Lwn values, although reduced, is still a problem in the 
coastal regions. More work is planned to improve the NIR 
methods and possibly will include absorbing aerosols in 
the aerosol model suite. Work is under way to improve the 
cloud detection and masking algorithm. Other algorithms 
will be examined and improved, such as the algorithm for 
the turbid water flag. 


28 



McClain, Barnes, Eplee, Jr., Franz, Hsu, Patt, Pietras, Robinson, Schieber, Schmidt, Wang, Bailey, and Werdell 


Chapter 4 


SeaWiFS Global Clear- Water Analysis 

Robert E. Eplee, Jr. 

SAIC General Sciences Corporation 
Beltsville, Maryland 

Charles R. McClain 
,V.4 5.4 Goddard Space Flight Center 
Greenbelt, Maryland 

Abstract 

The SeaWiFS CVT made a comparison of global dear-water radiances retrieved by SeaWiFS with normalized 
water-leaving radiances measured by MOBY as a check on the accuracy and stability of the vicarious calibration 
of SeaWiFS. The procedures and results of this comparison are described in this chapter. 


4.1 INTRODUCTION 

The SeaWiFS CVT performed a vicarious calibration 
of SeaWiFS (Eplee and McClain 2000a) using data from 
MOBY (Clark et al. 1997), deployed off of Lanai, Hawaii. 
As a check on the accuracy and stability of the vicarious 
calibration, the CVT made a comparison of global dear- 
water radiances, retrieved by SeaWiFS, with normalized 
water-leaving radiances (Lwn) measured by MOBY. For 
the purposes of this analysis, deep-water radiances are de- 
fined to be normalized water-leaving radiances (Gordon 
and Clark 1981) collected from ocean areas having a mini- 
mum depth of 1 km, while dear-water radiances are deep- 
water radiances collected from ocean areas with a max- 
imum chlorophyll a concentration of 0.15 mg m Radi- 
ances measured by MOBY met the clear- water criteria. 

4.2 CLEAR- WATER ANALYSIS 


CVT used the mean Lwn determined for MOBY during 
the vicarious calibration as the comparison for the Sea- 
WiFS clear- water radiances. 

The comparison of the two sets of radiances is shown in 
Table 8. The SeaWiFS dear-water radiances are slightly 
higher than the MOBY radiances because MOBY is lo- 
cated in water containing a marginally higher chlorophyll a 
concentration than that found for the global dear-water 
means. The SeaWiFS global dear-water radiances agree 
with the MOBY measurements of Lwn, verifying the va- 
lidity of the vicarious calibration of SeaWiFS. The appar- 
ent difference between SeaWiFS and MOBY in band 6 
(670 nm) occurs because radiance values near zero make 
the characterization of the vicarious gain difficult for this 
band (Eplee and McClain 2000a). For comparison, the 
nominal dear-water Lwn for the CZCS 520 and 550 nm 
bands are 0.495 and 0.280 mW cm 2 jim 1 sr respec- 
tively (Gordon and Clark 1981). 


The CVT computed global mean clear- water radiances 
for bands 1-6 from the time series of eight-day compos- 
ite SeaWiFS images produced by the Project using the 
calibration table and vicarious gains implemented for the 
third reprocessing. In this analysis, 124 eight-day compos- 
ites were used in the analysis, spanning a time range of 15 
September 1997 through 24 May 2000. 

As is discussed in Eplee and McClain (2000a), prob- 
lems with the surface ir radiance ( E s ) measurements for 
MOBY require the computation of Lwn f° r MOBY from 
measurements of water-leaving radiance {Lw), using the 
solar zenith angle and the atmospheric diffuse transmit- 
tance. Because the dear-water analysis requires the deter- 
mination of Lwn 3tt MOBY for dear-sky conditions, the 


Table 8. Clear- water radiance comparison between 
the MOBY vicarious radiances (in units of mW 
cm -2 sr -1 jzm -1 ) and the SeaWiFS dear-water ra- 
diances. The values shown are the arithmetic means. 


A 

[nm] 

MOBY 

Vicarious 

Radiance 

SeaWiFS 
Clear - Water 
Radiance 

412 

1.83536 

1.95934 

443 

1.61899 

1.75631 

490 

1.11291 

1.17869 

510 

0.651690 

0.697108 

555 

0.271239 

0.296487 

670 

0.0147519 

0.0348858 
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4.3 TIME-SERIES ANALYSIS 

The time series of the mean SeaWiFS clear-water and 
deep-water radiances are plotted in Figs. 12 and 13, respec- 
tively. The figures establish that the vicarious calibration 
and the time corrections for each of the SeaWiFS bands 
are stable over the course of the mission. The radiances for 
bands 1-3 in the two figures show the differences expected 
for clear-water and deep-water comparisons. Bands 1 and 
2 show slight periodicities, which may be a result of the in- 
terannual variability of global chlorophyll a concentrations 
(e.g., El Niho-La Nina) and due to the fact that the sam- 
pling of the ocean on a global basis shifts north and south 
seasonally. 


The stability of the vicarious calibration and time cor- 
rections is also demonstrated by the plot of the mean clear- 
water chlorophyll time series in Fig. 14. The periodic 
trends observed in the radiances for bands 1 and 2 are 
also apparent in the chlorophyll concentration. The mean 
deep-water chlorophyll time series plotted in Fig. 15 shows 
the effects of the El Niho-La Nina transition on the global 
chlorophyll abundance. 

The stability of the atmospheric correction over the 
course of the mission is shown by plots of the time series 
of the mean atmospheric correction parameter e(765,865), 
the mean Angstrom exponent at 510 nm, and the mean 
aerosol optical depth at 865 nm, shown for both the clear- 
water and deep-water analyses in Figs. 16-18, respectively. 
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Abstract 

The SeaWiFS CVT has looked for along-scan effects in SeaWiFS data by examining mean radiances of LAC 
scenes over MOBY and by analyzing the ratio of SeaWiFS-to-MOBY (S:M) match-up data. Analyses of the 
SeaWiFS along-scan data and S:M ratios show decreases in the water-leaving radiances retrieved by SeaWiFS as 
the scan angles or optical paths of the observations increase. These two analyses, which are independent of each 
other, both point to an overcorrection of the SeaWiFS data by the SeaWiFS atmospheric correction algorithm. 


5.1 INTRODUCTION 

SeaWiFS data can show variations in response as a 
function of pixel number within the scan. These along- 
scan effects have three primary sources: 

1) Variations in the instrument response as a function 
of pixel number within the scan (scan modulation 
effects), 

2) Variations in the atmospheric correction as a func- 
tion of optical path length, and 

3) Variations in ocean reflectance. 

The variations in instrument response were measured in 
the laboratory calibration of SeaWiFS (Barnes et al. 1994) 
and are accounted for by the scan modulation corrections 
applied during the calibration of the data. The analysis 
discussed in this chapter investigates possible variations in 
the atmospheric correction. 

5.2 ATMOSPHERIC EFFECTS 

The CVT analyzed variations in response as a func- 
tion of pixel number within the scan by averaging the scan 
lines in LAC scenes over MOBY on a pixel- by-pixel basis 
for all of the clear pixels in each of the scenes. The MOBY 
LAC data were selected to minimize areal variations in 
ocean reflectance while including coverage of a high qual- 
ity in situ time series. For this analysis, 149 scenes, span- 
ning the mission from 19 September 1997 until 4 February 
2000, were processed to level-2 and the resulting normal- 
ized water- leaving radiances were averaged for bands 1-6. 


The average chlorophyll a concentration for these data is 
approximately 0.1 mgm" 3 . The results of the analysis are 
mean radiances as functions of pixel number for the time 
period of the analysis. 

The mean Lwn values are plotted in Fig. 19. These 
plots show a possible trend with pixel number, especially 
for bands 1 and 2. The trend for band 1 is particularly 
evident in Fig. 20, where the Lwn values are plotted as 
a function of the secant of the scan angle at each pixel. 
The two sides of the scan (before and after nadir) are plot- 
ted separately. The GAC cutoff in the LAC data occurs 
at a scan angle of 45%; the first and last 146 pixels of 
LAC data are excluded from the GAC data sampling. Fig- 
ure 20 implies that the atmospheric correction of the data 
is not symmetric about nadir and that there is a rolloff 
in the measured radiances of approximately 10% at the 
GAC data cutoff. Examination of Figs. 19 and 20 suggests 
that the atmospheric correction algorithm overcorrects the 
radiances as the scan angle increases. 

To look for seasonal effects, the data were broken down 
into three periods of the year, based on the sequential day 
of the year. Period 1 covers days 1-99, and includes data 
from 1998, 1999, and 2000, and has 57 scenes. Period 2 
covers days 100-299, and includes data from 1997, 1998, 
and 1999, and has 54 scenes. Period 3 covers days 300-365, 
and includes data from 1997, 1998, and 1999, and has 38 
scenes. (Period 2 encompasses 200 days because sun glint 
during summer leads to a loss of data during this time.) 
The Lwn values in band 1 for these periods are plotted 
in Fig. 21. The different amplitudes of the radiances for 
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Fig. 19. The mean Lwn values as a function of pixel number along the scan line. 


Pre-Nadir 

Post-Nadir 



1.4 1-6 
Sec(Scan Angle) 


Fig. 20. The mean Lwn values as a function of the secant of the scan angle for each pixel 
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Fig. 21. The seasonal means of L W n values as a function of pixel number along the scan line. 
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Fig. 22. The S:M MOBY matchups for band 1 as a function of the secant of the scan angle. 
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Fig. 23. The S:M matchups for band 1 as functions of the secants of the solar and sensor zenith angles. 
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Fig. 24. The S:M matchups for bands 2-5 as functions of the secant of the scan angle. 
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each period are presumably due to seasonal variations in 
the chlorophyll concentration. The mean radiance for each 
period still shows a trend with pixel number. 

5-3 SCAN ANGLE OF MATCHUPS 

Because the along-scan analysis indicates an overcor- 
rection by the atmospheric correction algorithm as the 
scan angle increases, the CVT examined the SeaWiFS-to- 
MOBY (S:M) matchups discussed in Eplee and McClain 
(2000a) as functions of scan angle, solar zenith angle, and 
sensor zenith angle. The match-up ratios for band 1 are 
plotted as a function of the secant of the scan angle in 
Fig. 22. The band 1 ratios are plotted as functions of the 
secant of the solar and sensor zenith angles in Fig. 23. 
Similar effects are observed for bands 2-5. The match-up 
ratios for these bands are plotted as functions of the secant 
of the scan angle in Fig. 24. Each of these figures also show 
linear fits to the data. 

Examination of Figs. 22-24 shows a decrease in the 
SeaWiFS radiances as the various angles increase, or as 
the optical path through the atmosphere increases. The 
decrease is approximately 10% at scan angles and sensor 
zenith angles of 45%, which are correlated parameters. The 
magnitude of the decrease as a function of scan angle is 


comparable to the effect observed in the along-scan analy- 
sis. The decrease is approximately 7% for the solar zenith 
angle. The decrease in the match-up ratios as a function of 
the secant of the scan angle is comparable to the decrease 
in the water-leaving radiances observed in the along-scan 
analysis as a function of the secant of the scan angle. Be- 
cause the water-leaving radiances measured by MOBY are 
independent of the sensor geometry, the decrease in the 
match-up ratios with scan angle are due to decreases in 
the water-leaving radiances retrieved by SeaWiFS, as a 
function of scan angle. 

5.4 DISCUSSION 

Analyses of the SeaWiFS along-scan data and the S:M 
match-up ratios show decreases in the water-leaving radi- 
ances retrieved by SeaWiFS as the scan angles or optical 
paths of the observations increase. These two analyses, 
which are independent of each other, both point to an 
overcorrection of the SeaWiFS data by the SeaWiFS atmo- 
spheric correction algorithm. An alternative explanation, 
however, may be the bidirectional reflectance of the ocean 
surface, or f/Q effect (Morel and Gentili 1996), but pre- 
liminary evaluations were inconclusive. 
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Abstract 


In this chapter, a match-up procedure is described which compares the retrieved SeaWiFS aerosol optical 
thicknesses with data from in situ measurements. The aerosol optical thickness at 865 nm is a y-pro uc o 
the SeaWiFS atmospheric correction and is routinely retrieved from SeaWiFS measurements. This work is part 
of the SeaWiFS calibration and validation efforts in studying the aerosol optical properties over the ocean 
thereby validating aerosol models used in the atmospheric correction of ocean color sensors. The aerosol model 
is an integral part of the SeaWiFS atmospheric correction. The SeaWiFS aerosol retrieval algorithm, the data 
acquisitions from both SeaWiFS and the in situ measurements, and the match-up procedure are described. 
Finally, some preliminary comparison results are presented and discussed. 


6.1 INTRODUCTION 

Aerosols play an important role in climate forcing and 
biogeochemical cycling. They not only directly influence 
radiative transfer in the atmosphere and, hence, change the 
radiance reflected to space, but they also indirectly affect 
the radiation budget by providing cloud condensation nu- 
cleii that lead to cloud formation (Charlson et al. 1987 and 
1992). There have been continuous efforts in recent years 
with both ground in situ measurements and remote re- 
trieval of aerosol optical properties using aircraft and satel- 
lite sensors. The primary goals of the SeaWiFS mission 
are routine global ocean color measurements and ocean 
bio-optical property data (Hooker et al. 1992 and McClain 
et al. 1998). In retrieving the ocean near-surface signals 
from sensor-measured radiances at the satellite, however, 
the atmospheric effects must be removed. This is known 
as atmospheric correction , which removes more than 90% 
of the observed radiance in the visible spectrum (Gordon 
and Wang 1994a). 


The SeaWiFS atmospheric correction algorithm uses 
two NIR bands (765 and 865 nm) to estimate the aerosol 
optical properties and extrapolate these into the visible 
spectrum. Aerosol optical properties, in particular, aerosol 
optical thickness, is a by-product of the SeaWiFS atmo- 
spheric correction. Aerosol optical thickness at 865 nm, 
r a (865), is routinely retrieved from SeaWiFS measure- 
ments. The aerosol optical thickness is proportional to 
aerosol particle concentration, and is one of the most im- 
portant aerosol optical parameters. 

In this chapter, the effort to compare and validate the 
SeaWiFS aerosol optical products with in situ measure- 
ments, principally from Aerosol Robotic Network (AERO- 
NET) data (Holben et al. 1998), are outlined. Some other 
in situ measurements from calibration and validation cam- 
paigns within the SIMBIOS Project are also analyzed. 
There are two primary objectives behind these compar- 
isons. First, because r„(865) is part of the SeaWiFS stan- 
dard product suite, it warrants validation. With global 
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ground-based measurements of the r a data readily avail- 
able through AERONET. comparisons between the Sea- 
WiFS derived values and ground-based in situ measure- 
ments are very useful in evaluating the algorithm perfor- 
mance. The second objective of this work is to determine 
the validity of the suite of aerosol models currently used 
by the SeaWiFS Project for atmospheric correction. 

6.2 PROCEDURES 

In this section, the SeaWiFS aerosol optical thickness 
retrieval algorithm and its implementation in the SeaWiFS 
image processing system are described. Then, the algo- 
rithm is extended to retrieve the aerosol optical thicknesses 
in all of the SeaWiFS bands. Next, the data acquisition 
procedure and match-up criteria for both SeaWiFS and 
in situ observations are outlined. Finally, a data analysis 
strategy for both SeaWiFS and the in situ ground mea- 
surements is discussed. 

6.2.1 Theoretical Bases 

By defining the reflectance p = nL/(F 0 cos0 o ), where 
L is the radiance in the given viewing direction, F 0 is the 
extraterrestrial solar irradiance constant, and 6> 0 is the so- 
lar zenith angle. The total upward reflectance at the top 
of the ocean-atmosphere system, p t (A), measured at the 
two SeaWiFS NIR bands (765 and 865 nm) can be written 
as: 


(Gordon and Wang 1994a, and Wang and Gordon 1994), 
where p as (^j) is the single scattering aerosol reflectance at 
SeaWiFS band , which is given by 


Pas(^j ) 


7a (Aj) U)g(Xj) pa(0 , #0, A j) 
4 cos 0 cos Oq 


( 11 ) 


The parameters r 0 (Aj), u; 0 (Aj), and p a {0,9 o ,\j) are the 
aerosol optical thickness, the aerosol single scattering al- 
bedo, and the effective value of aerosol scattering phase 
function related to the single scattering case, respectively. 
The 6 term is the sensor zenith angles. 

In computing e(A a , A 7 ), \ 3 is usually taken at the longer 
NIR band, i.e., 865 nm for SeaWiFS. The value of c(A t , A j) 
characterizes the spectral variation of aerosol optical prop- 
erties. For a given solar and viewing geometry, e(A a ,A J ) 
depends only on the aerosol model; therefore, it forms the 
link between e(A a , \j) and the aerosol model. 

Using lookup tables, which were generated with a set 
of candidate aerosol models developed by Shettle and Fenn 
(1979), the p 0 (A) -I- p ra (A) values at the two SeaWiFS 
NIR bands can be converted to the single scattering re- 
flectance Pas(A), thereby providing e(765,865) values for 
given aerosol models. Hereafter, all e(765,865) values will 
be denoted as e with sub- and superscripts, as applicable. 
The SeaWiFS retrieved atmospheric correction parameter, 
e, was obtained by a weighted averaging over individual e 
values derived from a set of aerosol models. Two aerosol 
models with e~ and e + such that 


PtW = Pr{ A) + pa(A) -F p Ta (A) , (8) 

where p r (A), p a (A), and p ra (A) are the contributions from 
multiple scattering of air molecules (Rayleigh scattering), 
aerosols, and Rayleigh-aerosol interactions, respectively. 
Note that the surface sun glint and whitecap terms in the 
above equation have been ignored. This equation can be 
rewritten as 

PtW - Pt{ A) = pa( A) + Pra(A). (9) 

The value of the left hand side of (9) can be estimated 
from the sensor- measured radiance, p*(A), and the com- 
puted Rayleigh scattering reflectance, p r ( A), if the surface 
atmospheric pressure is known. This gives the aerosol and 
Rayleigh-aerosol interaction contributions of p a (A)+p ra (A) 
at the two SeaWiFS NIR wavelengths. By using a set of 
candidate aerosol models, the effects of the spectral varia- 
tion of the p a (A) T pro (A) values at the two NIR bands are 
then extrapolated into the SeaWiFS visible bands (Gor- 
don and Wang 1994a). The extrapolation was achieved 
through a process of aerosol model selection from evalu- 
ating the atmospheric correction parameter, e(A a , A ; ), de- 
fined as 

«<VA>) = (10) 

Pas\Aj ) 


r < e 


can be obtained, where e~ is for the aerosol model with 
the largest e value less than or equal to ?, and is for the 
aerosol model with the smallest e value greater than e. 

With the two retrieved aerosol models, the correspond- 
ing aerosol optical thicknesses (r 0 ) for a given wavelength 
A, r’(A) and r+(A), can then be estimated using (11). Fi- 
nally, the SeaWiFS aerosol optical thickness was obtained 
by interpolating between the two models as 

r a (A) = (l-r a )r a -(A) + r a r a + (A), (13) 

where 



is the interpolation ratio (r a ) between the two models. 

SeaWiFS routinely retrieves the aerosol optical thick- 
ness at 865 nm as a standard product. It is straight forward, 
however, to extend the current aerosol optical thickness re- 
trievals to the remaining SeaWiFS wavelengths using (13). 
After making a necessary interpolation for wavelengths 
that are slightly shifted from the SeaWiFS bands, compar- 
isons between SeaWiFS results and in situ measurements 
are possible. Both the Cimel sun-sky scanning radiome- 
ter and the MicroTops II sun photometer have spectral 
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Table 9. AERONET sites used for the aerosol match-up analyses. 


AERONET Station 

Latitude 

Longitude 

AERONET PI 

Bahrain 

26.32 

50.50 

C. McClain 

Bermuda 

32.37 

—64.70 

B. Holben 

Dry Tortugas 

24.60 

-82.80 

K. Voss and H. Gordon 

Kaashidhoo 

4.97 

73.47 

B. Holben 

Lanai 

20.83 

-156.99 

C. McClain 

San Nicolas Island 

33.26 

-119.49 

R. Frouin 


wavelengths at 440, 500, 670, and 870 nm corresponding to 
SeaWiFS bands 2 (443 nm), 4 (510 nm), 6 (670 nm), and 
8 (865 nm), respectively. The aerosol optical thicknesses, 
therefore, can be compared at these four wavelengths. 

It should be noted that the satellite and in situ data 
usually differ in both temporal and spatial characteristics. 
The in situ data are temporally averaged, while the satel- 
lite data are spatially averaged. 

6.2.2 SeaWiFS Data Acquisition 

The SeaWiFS aerosol optical thickness data were ob- 
tained by spatially co-locating a 25x25 pixel grid box 
around the pixel containing the ground-based measure- 
ment station, thereby providing a maximum of 625 Sea- 
WiFS retrievals in each matchup. A spatial homogene- 
ity (uniformity) test in the retrieved r a (865) value was 
then conducted to screen thin cirrus and high altitude 
aerosol contamination, because the in situ and satellite 
measurements are often looking through different atmo- 
spheric paths. Only those satellite data sets which passed 
the spatial homogeneity test were used for the match-up 
analyses. The SeaWiFS operational code has been modi- 
fied to output, at a pixel- by- pixel level, the aerosol optical 
thickness at a wavelength of 865 nm, values of retrieved 
two aerosol models, and the model partition ratio r a value. 
Aerosol optical thicknesses at all of the SeaWiFS wave- 
lengths can be calculated using (13). 

6.2.3 In Situ Data Acquisition 

The ground-based measurements, used for the aerosol 
optical thickness match-up analyses, come from two pri- 
mary sources: the automated Cimel sun-sky scanning ra- 
diometer managed as part of the AERONET network, and 
handheld MicroTops II sun photometers. In the following 
two sections, the data acquisition procedures from these 
two instruments are discussed. 

6.2.3. 1 In Situ Data from Cimel 

A select group of ground stations from AERONET was 
chosen. These instruments were located at either coastal or 
island stations and were operational for a reasonable length 


of time after SeaWiFS went into operation. Table 9 pro- 
vides the AERONET station name, location (latitude and 
longitude), and the corresponding responsible AERONET 
principal investigator (PI). Currently, efforts are underway 
to include additional AERONET stations. The retrieval 
of data from AERONET was automated to facilitate the 
match-up analyses. Once per month, a script is automat- 
ically run to access the AERONET database and retrieve 
r a data for the predetermined sites. 

For the match-up purpose, the ground-based measure- 
ments from AERONET were first reduced to include only 
those records within ±3h of the SeaWiFS overpass time 
for a given station. These records include the aerosol opti- 
cal thicknesses measured at the four spectral wavelengths 
(440, 500, 670, and 870 nm). As an initial quality control 
step, the data were averaged and some variation param- 
eters were computed to screen possible cloud contamina- 
tion. Only those data sets that have low temporal vari- 
ations (stable atmosphere) were then further reduced to 
±1 h of the SeaWiFS overpass time and used for the match- 
up analyses. Usually, the Cimel instruments routinely take 
one measurement every 15 min near local noon; therefore, 
for a given SeaWiFS file, there may be as many as eight 
AERONET measurements that qualify as a match for the 
2 h time window. Figure 25 shows the flow chart of the 
match-up procedure for the aerosol optical thickness ob- 
tained from SeaWiFS and Cimel measurements. The strat- 
egy in the validation study is not to compromise good data 
with bad data for the purpose of more matchups, i.e., it is 
preferred to screen out some good data to keep the high 
quality of data sets. 

6. 2. 3. 2 In Situ Data from MicroTops 

The hand-held MicroTops sun photometer data, were 
collected by various investigators in field campaigns asso- 
ciated with the SIMBIOS project. The cruises from which 
the MicroTops data were collected, along with the loca- 
tions and investigators, are listed in Table 10. Data from 
the MicroTops instruments are reprocessed from raw volt- 
ages using code adapted from the AERONET standard 
Cimel processing routines. This ensures that the r a data 
derived from the MicroTops measurements will be compa- 
rable to the data provided by AERONET. The number 
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Yes 



Fig. 25. A flow chart of the match-up procedure for the aerosol optical thickness obtained from SeaWiFS 
and Cimel measurements. 
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Table 10. MicroTops data used for the aerosol match-up analyses. 


Cruise 

Location 

Investigator 

MOBY Refit 1 

MOBY Station, Hawaii 

D. Clark 

Jason Project 2 

Sea of Cortez 

G. Feldman 

GoCal97-98 3 

Gulf of California 

J. Mueller 

JUL98NAN 4 

Massachusetts Bay 

B. Schieber and A. Subramaniam 

HOTS 5 

HOTS Station, Hawaii 

J. Porter and C. Motell 

HMS Rose 6 

East Coast of US 

G. Feldman 

USF-BNL 7 

TOTO, Bahamas 

M. Miller 


1) MOBY refit: Measurements obtained at the Marine Optical Buoy site at Hawaii. 

2) Jason Project: An expedition to the Sea of Cortez with Robert Ballard and the Jason Project. 

3) GoCal: A cruise to the Gulf of California (Sea of Cortez). 

4) JUL98NAN: A National Oceanic and Atmospheric Administration (NOAA)-sponsored cruise off Nantucket Island, Mas- 
sachusetts in July 1998. 

5) HOTS: Hawaiian Ocean Time Series cruise. 

6) HMS Rose : A cruise aboard Her Majesty’s Ship (HMS) Rose , a wooden tall sailing ship, which went from Miami, Florida 
to New York, New York. 

7) USF-BNL: A cruise to the Tongue of the Ocean (TOTO) by a group from the University of South Florida (USF) and 
Brookhaven National Laboratory (BNL). 
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Fig. 26. The retrieved SeaWiFS aerosol optical thicknesses r„(A) compared with the ground in situ mea- 
surements from the various AERONET stations for four wavelengths: a) 440 nm; b) 500 nm; c) 670 nm; and 
d) 870 nm. 
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Table 11. Five aerosol optical thickness, r 0 (865), samples of MicroTops data compared with SeaWiFS r a (865) 
data. r™ in , r™**, and r£ vg are the miniumum, manximum, and average from each set of measurements, respec- 
tively. 


Cruise 

_min 

'a 

MicroTops Data 

-max 

'a 

T avg 

'a 

SeaWiFS 

T avg 

' a 

GoCal97 

0.0196 

0.0821 

0.0342 

0.0632 

GoCal98 

0.0191 

0.0548 

0.0258 

0.0341 

Jason Project 

0.0285 

0.0727 

0.0433 

0.0297 

HMS Rose 

0.0638 

0.5829 

0.1500 

0.0692 

JUL98NAN 

0.0573 

0.1764 

0.0877 

0.0889 


of hand-held MicroTops measurements that match a given 
SeaWiFS file varies greatly because the measurement pro- 
tocol for these instruments is not well defined. In general, 
there should be a minimum of three MicroTops measure- 
ments per matched SeaWiFS file. The MicroTops data, 
however, tends to be much more variable than the Cimel 
data. The in situ data screening procedure, as outlined for 
the Cimel data set, was therefore not applied. 

6.3 PRELIMINARY RESULTS 

The SeaWiFS derived aerosol optical thicknesses with 
those from the ground in situ measurements were com- 
pared. As this is an ongoing research project, all results 
shown here are preliminary. Figures 26a-d provide an 
overall comparison results of r a (A) between SeaWiFS and 
Cimel measurements at the four wavelengths 440, 500, 
670, and 865 nm. The Cimel measurements were from 
the AERONET stations listed in Table 9. The number 
of data contributed to each plot in Fig. 6 from individual 
stations, from the top of the list to the bottom in Table 9, 
are 9, 8, 5, 2, 1, and 8. The station at Lanai, therefore, 
only contributed one point, whereas the Bahrain station 
provided nine points shown in Figs. 26a-d. The dotted 
lines in Fig. 26 are the 1:1 line. Although the compari- 
son results vary both in time and location, Fig. 26 shows 
that the comparisons agree reasonably well at the longer 
wavelengths (670 and 865 nm) for most of the AERONET 
stations. 

There is no obvious bias in the retrieved aerosol op- 
tical thicknesses. This implicitly indicates that the cali- 
bration at the SeaWiFS 865 nm wavelength is reasonably 
accurate. At the short wavelengths, however, it appears 


that SeaWiFS has the tendency of overestimating r a (A) 
with respect to the in situ measurements, in particular, for 
the relatively large aerosol optical thicknesses (r a > 0.15). 
This is most evident with the SeaWiFS r a (443) compar- 
ison results. Some possible sources that contributed to 
the comparison differences are from both satellite mea- 
surements (calibration, aerosol models, and cloud and thin 
cirrus contamination) and in situ data (mainly from instru- 
ment calibration and cloud contamination). More studies 
are needed to understand all of these. 

Similarly, the in situ MicroTops II data, which were 
from the various SIMBIOS calibration and validation cam- 
paigns, were compared with the SeaWiFS measurements. 
Table 11 shows five sample comparison results from five 
field experiments. For the MicroTops data in Table 11, 
the minimum r™ ln , the maximum r™ ax , and average r avg 
from each set of measurements are listed. There are high 
variations in the MicroTops measurements. Applying a 
data screen procedure for MicroTops data was not possi- 
ble, as was for the Cimel data, because none of MicroTops 
data were able to pass this data screening process. 
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Abstract 

Validation of SeaWiFS requires the use of in situ (field collected) data sets. The SeaWiFS and SIMBIOS Projects 
have sponsored numerous Pis to collect in situ optics and chlorophyll data for the purpose of comparing va ues 
to those derived from the SeaWiFS instrument. The match-up design described here uses field data stored in 
SeaBASS, match-up analysis software, and a plotting and statistics package to validate the primary SeaVViFS 
derived ocean products, i.e., normalized water-leaving radiance, chlorophyll a concentration, and the diffuse 
attenuation coefficient. 


7.1 INTRODUCTION 

A key method for validating SeaWiFS data is to com- 
pare the remote sensing values with coincident in situ mea- 
surements collected from ships or other oceanic platforms. 
NASA has sponsored field research activities to build a 
database of optical, pigment, and related in situ data for 
use in validating the derived products (McClain et al. 1992 
and 1998). These data, which are stored in SeaBASS, have 
been made available to those assisting in the SeaWiFS al- 
gorithm development and validation effort, e.g., the SIM- 
BIOS Science Team and the Moderate Resolution Imaging 
Spectroradiometer (MODIS) Oceans Team. 

Data from SeaBASS were employed in developing the 
first version of OC2 (OC2vl) during the first SeaWiFS 
Bio-optical Algorithm Mini-workshop (SeaBAM, O’Reilly 
et al. 1998) and the second version (OC2v2, O’Reilly et al. 
2000). The prelaunch strategy for the match-up compari- 
son process is described in McClain et al. (1996) and has 
been updated since launch (Hooker and McClain 2000). In 
this chapter, a brief overview of SeaBASS is followed by 


a discussion of the match-up procedures for the validation 
process. Finally, results of the current match-up analysis 
(October 2000) are presented. 

7.1.1 SeaBASS Data Archive 

The SeaWiFS Project built SeaBASS (Hooker et al. 
1994) as a local repository for in situ optical and pig- 
ment data to be used in satellite calibration and validation 
activities (e.g., in situ and satellite match-up analyses). 
The system is composed of two separate databases, a bio- 
optical archive with data from over 350 field campaigns, 
and a historical pigment archive which holds in excess 
of 300,000 fluorometric and high-performance liquid chro- 
motography (HPLC) pigment data records. Bio-optical 
data include worldwide measurements of apparent and in- 
herent optical properties from research vessels, moorings, 
drifters, and other platforms. Other oceanographic and 
atmospheric data, such as, estimates of primary produc- 
tivity and aerosol optical thickness, are also archived in 
SeaBASS. 
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Participants in the SIMBIOS Project are contractually 
obligated to submit data to SeaBASS, so the size of the 
archive has grown rapidly (McClain and Fargion 1999). 
SeaBASS also includes data contributions from the Sea- 
WiFS Field Program (Stanford Hooker, PI) and other US 
and foreign researchers. Access to SeaBASS is provided via 
the world wide web at http://seabass.gsfc.nasa.gov, 
and users may both browse the archive directories and 
query the bio-optical and pigment databases. To protect 
the publication rights of its contributors, access to Sea- 
BASS is limited (via a password) to members of the Sea- 
WiFS and SIMBIOS Science Teams and to other regularly 
contributing investigators who have been granted access 
on an individual basis. 

SeaBASS supports standard (two-dimensional) Amer- 
ican Standard Code for Information Interchange (ASCII) 
data files, which are easily managed from any computer 
platform and by most programming languages. This ap- 
proach also avoids the use of proprietary software on the 
part of a contributing investigator and allows easy integra- 
tion of the files into the world wide web. The architecture 
of a SeaBASS data file is simple: data are presented in 
columns (delimited by a space, tab, or comma) and pre- 
ceded by a series of predefined metadata headers. The 
headers provide descriptive information on the data file, 
such as date, time, location, column name and units, in- 
vestigators, and additional ancillary information. 

To assist with the standardization of SeaBASS data 
files, the Project developed feedback software and proto- 
cols to evaluate the format of submitted data files. A com- 
plete description of SeaBASS data formats, format veri- 
fication protocols, and quality control is provided at the 
SeaBASS web site at http://seabass.gsfc.nasa.gov. 

7.2 METHODS 

The procedure for comparing SeaWiFS image data to 
in situ data is described here, with particular attention 
given to the description of the exclusion criteria used. The 
procedure described, and the exclusion critera used, are 
designed to minimize the effect of bad data on the com- 
parison while not biasing the results. 

7.2.1 In Situ Match-up Data Files 

In situ data are compared to SeaWiFS image files by 
matching the two data sources in time and space. Each 
in situ data set is first summarized in a single ASCII file, 
known as a match-up file. The match-up file adheres to the 
aforementioned SeaBASS data format. The data included 
in the match-up file are L w ( A); E s ( A), either measured 
directly or extrapolated from profiled downwelling irradi- 
ance, Ed( 0 + )(A); C a \ and A"d(490). Location and mea- 
surement times are also necessary for each record in the 
match-up file. 

For all incoming data, proper descriptive documenta- 
tion and possibly, interactions with the data provider, are 


required to assure the usefulness of the in situ data for 
matchups. Many types of instrumentation and operational 
procedures produce data for SeaWiFS matchups. The ma- 
jority of the optical data come from well-known profiling 
radiometers, but a growing number of other approaches 
are being tested that collect data at a few discrete depths 
or from above the surface. 

7.2.2 The Match-up Procedure 

The SeaWiFS match-up procedure includes using both 
UNIX C-shell scripts (CSH) and Interactive Data Lan- 
guage (IDL) procedures. These are used to compare in- 
dividual SeaWiFS Hierarchical Data Format (HDF) files 
to in situ data records of individual measurements, i.e., 
stations. To generate the initial list of SeaWiFS files for 
consideration, the SeaWiFS Project HDF file database is 
queried for files that spatially and temporally match rec- 
ords in the in situ data source. The pixel and line number 
of the SeaWiFS file that matches the in situ location is de- 
termined. A region encompassing 101x101 pixels centered 
on the matched pixel is extracted from each of the result- 
ing initial level-la (LlA) files (GAC, LAC, and HRPT). 
These extracted data files are saved to a new file for local 
disk storage. 

For all valid LlA matches, a level -2 (L2) product is gen- 
erated from the extracted file for further analysis. A com- 
prehensive suite of SeaWiFS L2 products are generated, 
including all Lw{ A), Tuw(A), E s ( A), C a , A<*(490), and 
ancillary data (ozone, windspeed, and atmospheric pres- 
sure). In the generation of the L2 products, the Siegel NIR 
correction (Siegel et al. 2000) is applied, as well as the sun 
glint and out-of-band corrections. As with the extracted 
LlA files, the L2 files are saved to local disk storage. For 
each successful L2 file generated, a record including the 
path to the L2 file, the matched pixel and line number, 
and associated latitude and longitude are written to a file. 
This list file and the in situ match-up file are used as input 
to the IDL program where the match-up exclusion crite- 
ria are applied and the match-up plots and statistics are 
generated. 

7.2.3 Match-up Exclusion Criteria 

Only a small percentage of the candidate SeaWiFS files 
become final valid matches. A number of exclusion criteria 
have been formulated to provide an objective set of points 
for SeaWiFS validation which removes invalid or redun- 
dant data from consideration. The approach presented 
here is a result of numerous iterations. Nevertheless, fu- 
ture modifications to the existing set of exclusion logic is 
likely as a better understanding of both the SeaWiFS and 
in situ data are attained. 

The first exclusion criterion applied is a time difference 
between the in situ record and the satellite overpass. A 
time window of ±180 min from the satellite overpass is used 
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Local Time 


Fig. 27. Plots of hourly L w { 412) measurements from the BTBM over the course of a clear day: a) uncorrected 
Lwi. 412) measurements; and b) L\v (412) measurements corrected for the cosine of the solar zenith angle. 


because this window is the time period of reasonable illu- 
mination in most situations and, presumably, constant at- 
mospheric conditions. An analysis of Lw(412) normalized 
to the cosine of the solar zenith angle from the Bermuda 
Test-Bed Mooring (BTBM) throughout a clear day shows 
that a ±180 min window is reasonable (Fig. 27). 

In the current match-up set, the majority of candidate 
SeaWiFS files (i.e., those that coincide with in situ sta- 
tions), are eliminated from consideration as a result of the 
SeaWiFS image pixels being flagged (excluded), usually 
because of clouds and stray light, although other factors 
are often present. Pixels are excluded if any of the fol- 
lowing flags are applied: atmospheric correction failure, 
land, sun glint, total radiance above the knee value, high 
satellite zenith angle, stray light, clouds or ice, coccol- 
ithophores, and low Lvw(555). Currently, pixels are not 
excluded with the turbid water flag applied. The authors 
have found that the current turbid water flag is not robust. 
It may erroneously flag pixels near cloud edges and fail to 


set in extremely turbid water if the SeaWiFS derived chlo- 
rophyll concentration is high (greater than approximately 
lOmgm -3 ). 

If a matchup passes the temporal exclusion criterion, 
the valid SeaWiFS pixels are averaged for a region en- 
compassing the matched pixel, and simple statistics are 
recorded. The match-up approach uses a 3x3 pixel box 
for this spatial comparison. It should be noted that Sea- 
WiFS GAC files have a lower sampling resolution than the 
LAC or HRPT products (approximately 4 km samples ver- 
sus approximately 1km samples, respectively). This sam- 
pling aspect results in a poorer representation of in situ 
conditions. 

Once the match-up code has reported all unflagged 
matches, the routine reviews the records for further re- 
finement. Calculations currently used for the exclusion of 
points and the match-up values include: 

1. Minimum number of valid pixels: At least five out of 

nine pixels in the 3x3 box considered from the Sea- 
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WiFS image must be valid (unflagged) to avoid con- 
taminated data. 

2. SeaWiFS file reduction: If there are GAC, onboard 
LAC, and HRPT matches corresponding to the same 
in situ point, the order of preference is onboard LAC, 
HRPT, and GAC. When multiple files of the same pre- 
ferred type match a single in situ record, the closest 
temporal match is selected. This selection is neces- 
sary when multiple HRPT stations or multiple GAC 
swaths cover the in situ point. Permanent real-time 
HRPT stations are given precedence over other HRPT 
sources. 


Values of E s are usually measured directly (above the sur- 
face), but are often estimated by extrapolating subsurface 
measurements of Ed through the air-sea interface, denoted 
as Ed( 0 + , although, there may be errors introduced in do- 
ing so (D. Clark, pers. comm.). The second calculation of 
Lwn (item c in the list above) is determined using the Lw 
values and a theoretical estimate of E a (M. Wang, pers. 
comm.): 


E S (X) = F o (0 o )cos(0 o )exp 


7 

cos(0 o ). * 


where 


(15) 


3. Duplicate in situ data reduction: Cases where multiple 
in situ casts are performed at the same station are re- 
duced to one representative record before the matchup 
by selecting the cast with the highest Lw(490) value 
normalized to the theoretical E s value. 

4. Multiple in situ measurements for each SeaWiFS file: 
Along-track measurements are included in the in situ 
match-up files where adjoining measurements are grea- 
er than 12 km apart. 

5. Surface irradiance calculation: Theoretical E s values 
are calculated using time and location for each in situ 
data record. This value can be used for later Lwn 
analysis. 

6. Solar elevation correction: A cos(#o) correction is ap- 
plied to L w match-up points (in situ and satellite). 
This adjustment helps remove time dependencies due 
to changes in the solar elevation and allows for a larger 
time window to be used. 

7. Out-of-range /Q( 4 90) values: Cases where the /Q(490) 
values are below the /Q(490) value of pure water (0.016 
in -1 ) or above a value of 6.4 m -1 are removed. 

8. Large coefficient of variation elimination: Satellite 

matchups with extreme variation between pixels in the 
SeaWiFS 3x3 pixel box (coefficient of variation, or the 
standard deviation divided by the mean value, greater 
than 0.2) are excluded. These typically represent fron- 
tal regions or other anomalies (e.g., cloud edges) in the 
SeaWiFS imagery, which make the match-up validity 
questionable. The coefficient of variation value was se- 
lected because the satellite data appeared to form sep- 
arate clusters of both fairly uniform data and highly 
variable data with a gap between them at about 0.2. 
Presently, this test is applied to each geophysical field 
(Lwn, C a , etc.) independently. 

Currently, three comparisons can be performed for each 

candidate matchup: 

a) SeaWiFS Lw to in situ Lw\ 

b) SeaWiFS Lwn to in situ Lwn using observed Lw 
and E s (Lwn = F$Lw /E s \ and 

c) SeaWiFS Lwn to in situ Lwn using observed Lw , 
and theoretical E s . 


1 = + fitcMA) + 0.0054, (16) 

A is the wavelength, 6 0 is the solar zenith angle, r r is the 
Rayleigh optical thickness, a oz is the ozone absorption co- 
efficient, and ft tc is the total columnar ozone in Dobson 
units (DU). The wavelength-independent constant value 
of 0.0054 accounts for the effects of aerosols and was es- 
timated using an aerosol optical thickness of 0.1 and the 
M90 aerosol model. The ozone value used is an average 
of the ozone values for the 3x3 pixel box surrounding the 
matched pixel. This same ozone value is the one used in 
the operational processing. 

An extensive comparison of both in situ and theoretical 
clear-sky E s values for the valid matchups was conducted. 
The comparison showed reasonably good agreement be- 
tween the measured and theoretical £,(A) values. There 
were, however, several values above the maximum theoret- 
ical value. These may be the result of cloud reflection ef- 
fects, poor calibrations, etc. The theoretical E s value was 
shown to give similar match-up comparisons than those 
using the in situ E s values. The theoretical E a values, 
therefore, are used to compute the Lwn values for the in 
situ data. 

Previous versions of the exclusion criteria eliminated 
match-up points in cases where the measured (in situ) E a 
values differed substantially from the theoretical value in 
(15). This step was originally included to remove matches 
that may have been contaminated by clouds, but were not 
flagged as such in the level -2 processing. This step was 
later removed as it was noted the approach does not gener- 
ally exclude any additional points from the final matchups. 
Other steps in the exclusion processing removed faulty 
records. Another test was to compare the difference in so- 
lar zenith angles between the satellite and in situ measure- 
ment times for Lw measurements. This condition, how- 
ever, was ameliorated by applying the cos(0o) correction 
for Lw comparisons as mentioned above, so this exclusion 
criteria was also dropped. 

Once all the exclusion criteria are applied, statistics 
are generated for each Lw( A), Lwn( A), C a , and ZQ(490) 
comparison. These are written to an ASCII file. The final 
matched data set is saved as an IDL saveset (Table 12) for 
future reference and analysis. 


48 










SeaWiFS Postlaunch Calibration and Validation Analyses, Part 2 



In Situ Chlorophyll a 



Fig. 28. (cont.) Comparisons of SeaWiFS and in situ measurements of Lwn and f) chlorophyll a, and g) 
N(490). Negative Lwn values have been excluded. 


Table 12. Fields stored in a match-up results summary file. 


Name 

Description 

CRUISEID 

Cruise identifier 

ENV.YEAR 

Year of in situ measurement 

ENV .MONTH 

Month of in situ measurement 

ENV-DAY 

Day of in situ measurement 

ENV.H0UR 

Hour of in situ measurement 

ENV -MINUTE 

Minute of in situ measurement 

ENV.SECOND 

Second of in situ measurement 

ENV.LAT 

Latitude of in situ measurement 

ENV-LON 

Longitude of in situ measurement 

ENV_LW(A) 

In situ Lw at wavelength A 

env_es(A) 

In situ E a at wavelength A 

ENV.KD490 

In situ Kd( 490) 

ENV.CHL 

In situ chlorophyll concentration 

TDIFF 

Time difference between in situ measurement and satellite overpass 

SAT-FILE 

Satellite extract file name 

SAT-COUNTS 

Number of valid pixels in 3x3 box surrounding in situ lat/lon 

SAT-SOLZ 

Solar zenith angle at satellite overpass 

SATJSOLA 

Solar azimuth angle at satellite overpass 

SAT-SENZ 

Satellite zenith angle at satellite overpass 

SAT-SENA 

Satellite azimuth angle at satellite overpass 

SAT-LW-(A) 

Satellite Lw at wavelength A 

SAT-ES-(A) 

Satellite E s at wavelength A 

SAT-NLW_(A) 

Satellite Lwn at wavelength A 

SAT-TAUA-(A) 

Satellite aerosol optical thickness at wavelength A 

SAT_ANGSTR0M_555 

Satellite Angstrom exponent between 865-555 nm 

SAT.EPSILON 

Satellite e value 

SAT-OZONE 

Ancillary ozone value 

SAT-WINDSPEED 

Ancillary wind speed 

SAT-PRESSURE 

Ancillary atmospheric pressure 

SAT-CHL.0C2 

Satellite-derived chlorophyll using OC2v2 (operational algoithm at the 
time of this writing) 
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Table 12. (cont.) Fields stored in a match-up results summary file. 


Name 

Description 

SAT_K_490 

Satellite-derived ^(490) 

SAT_(A)_STATS 

Simple statistics on satellite measurements (minimum, maximum, 


and standard deviation) 

ENV_ES(A)_CALC 

Calculated in situ E s using the Gordon clear sky model (Gordon et al. 1988) 


Table 3. Decrease in the number of matches as different criteria are applied, in order, to the match-up 


set. 

Processing Step 

Number of Matches Remaining 


Initial set (at least one valid SeaWiFS pixel) 

886 


Five out of nine pixels valid 

659 


Duplicate satellite coverage 

361 


Duplicate in situ measurements 
Coefficient of determination tests 
(negative Lw values excluded): 

198 


Lw(412) 

119 


Lw(443) 

141 


Lw(490) 

153 


Lw(510) 

136 


Lw (555) 

149 


C a 

76 


K d { 490) 

125 


Table 14. Lw statistics are presented for the first five wavelengths and C a values. The mean ratios are 
the SeaWiFS- to-in situ values, a is the standard deviation, and CV is the coefficient of variance. The 
number of points for each parameter is not the same because of differences in the specific wavelengths 
measured in situ by different investigators. C a values were sometimes measured without coincident 
optics, or, in other cases, the reverse was true. Some matchups were dropped because of negative L w 


values. 


vouuco. 

Parameter 

Mean 

Patio 

a 

CV 

Number of 
Observations 

Range of In Situ 
Values 

L w ( 412) 

0.856 

0.222 

0.259 

119 

0.330-3.238 

L w ( 443) 

0.957 

0.246 

0.256 

141 

0.219-2.709 

Lw(490) 

0.973 

0.206 

0.211 

153 

0.252-1.695 

Lw(510) 

1.075 

0.236 

0.219 

136 

0.215-0.857 

Lw(555) 

1.106 

0.279 

0.252 

149 

0.127-0.538 

c„ 

1.061 

0.498 

0.434 

76 

0.062-4.650 

K d ( 490) 

1.255 

0.287 

0.222 

125 

0.018-0.340 


7.2.4 Match-up Results 

The current (October 2000) in situ data records in Sea- 
BASS were processed for each exclusion criterion outlined 
above (Table 13). Of the 2,389 in situ records considered 
(39 data sets and 4,248 SeaWiFS files), 198, or 8.3% passed 
all exclusions to become final matches. 

Comparisons of Lwn between in situ sites and Sea- 
WiFS show generally good agreement in the radiance val- 
ues (Figs. 28a-e and Table 14), especially at the higher ra- 


diance values. Note the lower number of match-up points 
in the Lwn{ 412) and, to a lesser extent, Lwn( 443) band. 
This reduction at lower wavelengths occurs in coastal ar- 
eas where the spectral water-leaving radiance in the NIR 
wavelengths is greater than zero, a condition that causes 
the overestimation of aerosol radiances and low, often neg- 
ative, Lw values in these bands. 

Negative Lw values can also result when absorbing 
aerosols are present with the most severe effects being at 
the shortest wavelengths. In nearly all cases, this anomaly 
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occurred in Case-2 waters or in a phytoplankton plume 
with relatively high C a concentrations. Many of these val- 
ues exhibited high coefficients of determination (greater 
than 0.5) in the 3x3 pixel set. Note that the lower number 
of points in the Lw{ 510) panel is not due to overestima- 
tion of the aerosol radiance, but stems from the absence of 
this channel in some field radiometers. 

The comparison of Lwn(412) (Fig. 28a) shows that 
SeaWiFS retrievals have a slight negative bias, which in- 
creases at lower values. The conditions reported above 
that resulted in negative Lwn retrievals are influencing 
the retrievals at low values as well. The Lwn (443) com- 
parison (Fig. 28b) does not show the same bias seen in 
the Lwn (41 2) comparison, although it is possible that the 
effect is present at this band, except at a much smaller 
extent. SeaWiFS compares quite well with the in situ 
measurements for Lwn (490), Lwn (510), and Lwn(555) 
(Figs. 28c-e, respectively), although Lwn (555) may be 
slightly positively biased. 

The comparison of SeaWiFS and in situ chlorophyll 
values (Fig. 28f) demonstrates the reasonably high quality 
of the chlorophyll a data product for m situ values from 
about 0.02-7 mg m~ 3 . For this comparison, SeaWiFS de- 
rived chlorophyll a values were determined using a different 
ocean chlorophyll algorithm (OC4v4). A comparison of the 
two algorithms showed that the OC4 algorithm performed 
better than the OC2 algorithm. OC4 became the opera- 
tional algorithm with the third reprocessing of SeaWiFS 
data. These results are affected not only by variations 
in SeaWiFS imagery, but from inaccuracies in field chloro- 
phyll measurement techniques. Ongoing work continues to 
improve the SeaWiFS and in situ comparisons from both 
the SeaWiFS imagery and field measurement aspects of 
the effort. 

7.3 CONCLUSIONS 

The SeaWiFS validation effort has demonstrated the 
capability of matching SeaWiFS and in situ measurements. 
The comparison of SeaWiFS and m situ chlorophyll con- 
centration demonstrates that the SeaWiFS C a values are 
reasonably good. A weakness in the validation effort has 
been the quality of the in situ data used as truth against 
the SeaWiFS products. Problems of instrument deploy- 


ment, calibration, pigment estimation, data processing, 
and documentation have resulted in errors leading to in- 
valid matchups despite the calibration and data analysis 
round-robins, protocol development, quality control, and 
other activities initiated by the SeaWiFS Project to im- 
prove data quality. Improvements, however, are being re- 
alized as stricter data ingest methods (i.e., FCHECK) and 
other activities proceed. 

The SeaWiFS and in situ match-up effort has provided 
a quantitative evaluation of the data products, although, 
it is limited to regions where data have been reported. At 
present, the global representation of Case-1 (open ocean) 
waters is not complete, especially in areas where chloro- 
phyll a concentrations are less than 0.05 mg m~ 3 and in 
polar regions. Figure 29 provides the locations of the cur- 
rent match-up data set, which corresponds to the points in 
Fig. 28. The SeaWiFS and SIMBIOS Projects are actively 
supporting continued field research in these undersampled 
areas. 

Efforts to report individual match-up results from the 
ongoing validation effort on-line have been undertaken to 
allow individual investigators the opportunity to assess 
their data further and perhaps explore methods for future 
in situ data improvements. It is hoped this exchange of 
ideas will continue to improve the accuracy and usefulness 
of the validation data set. In addition, other approaches to 
validating pigment concentrations are being explored, e.g., 
comparing regional and seasonal statistics where sufficient 
in situ data are available. As the ongoing calibration and 
validation efforts expand under SIMBIOS and other mis- 
sions, and as improvements in the atmospheric correction 
and bio-optical algorithms are developed, the accuracy of 
the products will continue to improve. 

Future enhancements planned for the match-up analy- 
sis include the addition of an exclusion criterion based on 
the optical depth of the water column. If the optical depth 
is greater than the depth of the water column, bottom 
reflection will influence the Lwn values retrieved by Sea- 
WiFS. Currently, there is no way of excluding these points 
from the match-up data set. Plans are also underway to 
improve the SeaWiFS turbid water flag. Should a suffi- 
ciently robust flag be developed, the match-up data set 
will be partitioned into Case-1 and Case-2 water types. 
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Glossary 

AERONET Aerosol Robotic Network 
AI Absorbing Aerosol Index 
AOT Aerosol Optical Thickness 
ASCII American Standard Code for Information In- 
terchange 

BNL Brookhaven National Laboratory 
BTBM Bermuda Test Bed Mooring 

CSH UNIX u C-shell” script programming utility 
CVT Calibration and Validation Team 
CZCS Coastal Zone Color Scanner 

DU Dobson Unit of total ozone 

GAC Global Area Coverage 
GoCal Gulf of California 

HDF Hierarchical Data Format 
HMS Her Majesty’s Ship 
HPLC High-Performance Liquid Chromotography 
HRPT High Resolution Picture Transmission 

IDL Interactive Data Language 

JUL98NAN A NOAA-sponsored cruise off Nantucket Is- 
land, Massachusetts in July 1998. 

LI Level -1 SeaWiFS data product 
L1A Level-la SeaWiFS data product with naviga- 
tion information 

L2 Level -2 SeaWiFS data product 
L3 Level-3 SeaWiFS data product 
LAC Local Area Coverage 

MOBY Marine Optical Buoy 

MODIS Moderate Resolution Imaging Spectroradio- 
meter 

NASA National Aeronautics and Space Administra- 
tion 

NIR Near-Infrared 

NOAA National Oceanographic and Atmospheric Ad- 
ministration 

OC2 Ocean Chlorophyll 2 algorithm 
OC2vl OC2 version 1 
OC2v2 OC2 version 2 

OC4 Ocean Chlorophyll 4 algorithm 
OC4v3 OC4 version 3 

PI Principal Investigator 

SeaBAM SeaWiFS Bio-optical Algorithm Mini- 
workshop 

SeaBASS SeaWiFS Bio-optical Archive and Storage Sys- 


Symbols 

a 0 * Ozone absorption coefficient. 

b Intercept. 

C a Chlorophyll a. 

C a Mean chlorophyll a values. 

Cp Phaeophytin a. 

CV Coefficient of variance. 

Ed( A) Downwelling irradiance. 

Ed(0 + ) Extrapolated subsurface measures of Ed{ A) through 
the air-sea interface. 

F a (A) Surface irradiance incident at the ocean surface, 
E d { 0 + ). 

/ Scan angle. 

Fo(A) Solar irradiance constant adjusted for the Earth- 
sun distance variations. 

f/Q Ocean surface bidirectional reflectance parameter. 

i Square root of — 1 . 

I Aerosol index. 

K (490) Diffuse attenuation coefficient of seawater measured 
at 490 nm. 

A'd(490) Diffuse attenuation coefficient of downwelling irra- 
diance at 490 nm. 

L(A) Radiance. 

L a (A) Radiance measured at the TOA from aerosol scat- 
tering alone. 

(L a (Aj)) Aerosol interaction term. 

L r (A) Radiance measured at the TOA from Rayleigh scat- 
tering alone. 

(L r (Ai)) The Rayleigh contribution. 

L ra (A) Radiance at the TOA from Rayleigh- aerosol inter- 
active scattering. 

(L ra (Aj)) Rayleigh- aerosol interaction term. 

Lt(A) Radiance measured at the TOA. 

(L t (\i)) The radiance at the TOA measured by the SeaWiFS 
instrument. 

Lw(A) Water-leaving radiance. 

Lwn{ A) Normalized water-leaving radiance. 

Lwn( A,) Normalized water-leaving radiance for a given band. 

L'wn( A«) Normalized water-leaving radiance at the band- 
centered wavelength. 

(L(Aj)) Average radiance for band i as weighted by the sen- 
sor spectral response function. 

{{LwN(Ai))} Average of the parameter (L(Aj)) as weighted by 
the sensor spectral response function. 

m Slope. 

pa (A) The effective value of aerosol scattering phase func- 
tion related to the single scattering case. 


tern 

SeaWiFS Sea-viewing Wide Field-of-view Sensor 
SIMBIOS Sensor Intercomparison and Merger for Bio- 
logical and Interdisciplinary Oceanic Studies 
SRF Spectral Response Function 

TOA Top of the Atmosphere 
TOMS Total Ozone Mapping Spectrometer 
TOTO Tongue of the Ocean study (Bahamas) 

USF University of South Florida 


R Correlation coefficient. 

R\\l The TOA reflectance ratio measured by SeaWiFS. 
Rill The reflectance ratio calculated assuming an atmo- 
sphere containing only large nonabsorbing particles 
and Rayleigh scattering molecules. 
r(A*) Ratio value of Lwn at the SeaWiFS band-centered 
wavelength versus the average values weighted by 
the SeaWiFS spectral response function. 
r 0 (p) Mode radius of the particle size distribution in mi- 
crons. 
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r a Interpolation ratio between two models in the de- 
rived epsilon values. 

ri(Xi) Ratio value of Lwn at the SeaWiFS band-centered 
wavelength versus the average values weighted by 
the SeaWiFS spectral response function. 

R n ( A) Remote sensing reflectance. 

Rtb(X) Measured remote sensing reflectance. 

Rm(X) Calculated remote sensing reflectance. 

R^ A compact notation to represent R rs (X i ) / R rs (Xj) in 
the definition for aerosol absorbing index, see (1). 

S*(A) The SeaWiFS spectral response function for band i. 

i(A) Diffuse transmittance of the atmosphere at the sen- 
sor viewing direction. 

t(Aj) Atmospheric diffuse transmittance at the instrument 
viewing direction. 

<o(A) Diffuse transmit tance of the atmosphere at the solar 
direction. 

x The abscissa. 

y The ordinate. 

7 See (16). 

c Ratio of the single-scattering aerosol radiances. 

? Weighted average of individual e values. 

Aerosol model with the largest e value greater 
than or equal to ?. 

€' Aerosol model with the smallest t value greater 
than e. 

e(Ai, Xj) Atmospheric correction parameter defined as the ra- 
tio of single-scattered aerosol reflectance between 
two wavelengths, p as (Xi)/p a 9 (Xj). 

0 Viewing zenith angle measured at the pixel. 

0o Solar zenith angle measured at the pixel. 

A Wavelength. 

A i, Xj Wavelengths corresponding to sensor spectral bands 
i and j. 

p Reflectance. 

p 0 (A) Reflectance measured at the TOA from aerosol scat- 
tering alone. 

paa (A) Single-scattered aerosol reflectance. 

p r (A) Reflectance measured at the TOA from Rayleigh 
scattering alone. 

Pro (A) Reflectance at the TOA from Rayleigh- aerosol in- 
teractive scattering. 

p t (A) Reflectance measured at the TOA. 

Pw(X) Water- leaving reflectance at the sea surface. 

a Standard deviation. 

o' Width of the particle size distribution. 

r a (A) Aerosol optical thickness. 

r* vg Average r a ( A) from each set of measurements. 

T max Maximum r 0 ( A) from each set of measurements. 

r™ ,n Minimum r 0 (A) from each set of measurements. 

r r Rayleigh optical thickness. 

u> 0 (A) Aerosol single scattering albedo. 

Dtc Total columnar ozone in Dobson units (DU). 
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